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Abstract 

Liischer's method is routinely used to determine meson-meson, meson-baryon and baryon-baryon 
s-wave scattering amplitudes below inelastic thresholds from Lattice QCD calculations - presently 
at unphysical light-quark masses. In this work we review the formalism and develop the requi- 
site expressions to extract phase-shifts describing meson-meson scattering in partial-waves with 
angular-momentum I < 6 and / = 9. The implications of the underlying cubic symmetry, and 
strategies for extracting the phase-shifts from Lattice QCD calculations, are presented, along with 
a discussion of the signal-to-noise problem that afflicts the higher partial-waves. 
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I. INTRODUCTION 



The s-wave interactions between hadrons are being calculated with Lattice QCD (LQCD) 
with increasing precision. Presently, such calculations are being performed at unphysical 
light-quark masses, and in the case of mesonic interactions, extrapolations to the physical 
light-quark masses are made possible by chiral perturbation theory (xPT). Unfortunately, 
such extrapolations are presently not reliable for baryon-baryon interactions, and it is likely 
that LQCD calculations at, or very near, the physical light-quark masses will be required 
to make precise predictions for these interactions due to the fine-tunings that are known to 
exist in nuclear physics. In most LQCD calculations, periodic boundary-conditions (BC's) 
are imposed on the quark and gluons fields in the spatial-directions of the lattice volume 
and Liischer's method ^ can be used to extract scattering phase-shifts from the energy- 
eigenvalues of two-hadron states that lie below inelastic thresholds. As it is the irreducible 
representations (irreps) of the cubic group that determine the degeneracies of the eigenstates 
in the (cubic) lattice volume, it is difficult to determine the phase-shifts, 5i, beyond the 
lowest few partial-waves. Each of the irreps of the cubic group have a non-zero overlap with 
infinitely many irreps of S0(3), and as a result, the energy-eigenvalues of two-hadron states 
transforming as a certain irrep of the cubic group receive contributions from the phase- 
shifts in an infinite number of partial- waves. In contrast, two-particle systems confined 
in a harmonic oscillator potential have a one-to-one relation between phase shifts and the 
energy-eigenvalues since the potential respects SO (3) symmetry [SI S] ^- The mixing of 
angular momentum in cubic irreps consequently limits the precision with which the phase- 
shift in any given partial-wave can be extracted in a LQCD calculation. This was made 
obvious in the work of Mandula, Zweig and Govaerts,[6J, and explicitly detailed in Liischer's 
papers [H E]. Liischer calculated the energies of states in the irrep of the cubic group as 
a function of 5o and ^4, and gave general expressions for the energies of states transforming 
in each of the cubic irreps in terms of the 5i. The extension of this formalism to systems with 
non-zero center-of-mass (CM) momentum was performed by Rummukainen and Gottlieb [7], 
and later by Kim, Sharpe and Sachrajda [8]. Recently finite- volume expressions for three- 
nucleon systems within cubic volumes have been investigated [HHTT]. 

Fully-dynamical Uf = 2+1 LQCD calculations of meson-meson interactions in the isospin- 
stretch-states {i.e. no disconnected diagrams) are presently enabling predictions of the s- 
wave interactions with percent-level precision [T2HT7] (for a recent review, see Ref. [T8]). 
and very recently a preliminary calculation of the Tr+Tr"*" d-wave phase-shift has been per- 
formed [20]. Further, preliminary calculations of / = vrvr scattering, which contain discon- 
nected diagrams, have been performed [21]. These calculations were preceded by quenched 
LQCD calculations [22| - H0] . and by early rif = 2 LQCD calculations [H]. Meson-baryon 
systems are starting to be explored in the channels for which disconnected diagrams are 
not required in the LQCD calculations [12]. Further, LQCD calculations of baryon-baryon 
interactions are beginning to become reliable at unphysical pion masses [13H1S], and recently 
the binding energy of the H-dibaryon has been calculated [17]. Now that the methodology 
for extracting s-wave interactions has been shown to be effective, it is appropriate to explore 
the higher partial- waves. In the meson sector, a determination of the p-wave phase-shifts 
has direct implications for post-dieting the mass and width of the p- meson [HI |19] , but this 



The same is true when a "spherical-wall" is imposed on the separation between hadrons, as has been 
demonstrated in recent lattice effective field theory (LEFT) calculations [5]. 
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requires evaluating disconnected diagrams in LQCD - calculations that are computationally 
expensive. In nucleon-nuclei scattering, the experimentally determined p-wave phase-shifts 
are thought to be at the heart of the "A^-puzzle" in nucleon-deuteron scattering. Further, 
it is found phenomenologically that only the phase-shifts in partial waves with 4 are 
required to perform relatively precise calculations of nuclear structure and reactions (at the 
physical pion mass). 

The formalism required to analyze the J = 1 coupled-channels, in which the deuteron is 
the ground-state, has been put in place by Liu, Feng and He [52], and exploratory quenched 
calculations of the s-d mixing parameter, e, at pion masses of m.„ ~ 730,530 and 380 MeV 
have been performed in a small number of lattice volumes Further, there has been 

recent work in developing the phenomenology that goes beyond Liischers formalism in an 
attempt to explore resonances (and couplings to multi-hadron final-states) in the single 
baryon and meson sectors [52]. 

It is appropriate to point out that there is a substantial amount of information and 
technology that is directly relevant to this subject, in particular space-groups, that has been 
developed for study of condensed matter systems. Much of the work in this paper draws 
directly from various applications found in these fields. Discussions of space-groups can be 
found in texts, such as Ref. [S3] or Ref. [51], as are discussions of point-groups, and other 
formalisms that impact the present calculations. 

While the papers by Liischer [H |2] contain the required formalism, we take this opportu- 
nity to present the explicit relations between the energy-eigenvalues of two-meson states in 
a cubic volume and the phase-shifts in the partial- waves with / < 6 and 1 = 9. The exper- 
imentally measured phase-shifts describing tttt scattering in the lowest-lying partial-waves, 
appropriately parameterized, are used to perform estimates of the energy-eigenvalues that 
are expected in LQCD calculations of such systems over a range of lattice volumes. We 
also discuss the issue of signal-to-noise degradation while performing lattice calculations in 
higher partial waves. 

II. FORMALISM 

In the absence of interactions, the states in the cubic volume can be defined by their behavior 
under transformations of the cubic group and by their energy. As the momentum in the 
volume is quantized in integer multiples of 27r/L, where L is the spatial-extent of the volume, 
the energy quantum number can be replaced by the magnitude of the integer triplet defining 
the momentum, |np, where n = {nx,ny,nz). Instead of the energy, it is convenient to refer 
to the particular |np-shell. For each partial-wave with / < 6 and / = 9, an irrep of the cubic 
group is identified for which 6i provides the dominant contribution to the interaction energy. 
Sources and sinks used in LQCD calculations that are constructed to transform under such 
irreps will allow for a determination of the 6i at some level of precision. The energy of 
states with |np < 6 are required to lie below the inelastic threshold in order to obtain all of 
the phase-shifts with / < 6, thereby requiring relatively large lattice volumes. Further, the 
energy of a state in the |np = 14-shell is required to obtain the I = 9 phase-shift. 

A non-zero phase-shift in a given partial-wave will, in general, contribute to the energy- 
eigenvalues of two-hadron states in the volume that transform as one or more irreducible 
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TABLE I: Decomposition of the orbital angular momentum eigenstates, \l,m), into irreps of the 
cubic group, F^*), for Z < 9 (see, for instance, Ref. [53j). 
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representations of the full cubic group, T^'^\ ^ Table [l] shows the decomposition of the orbital 
angular momentum eigenstates, \l, m), into the T^*^ for / < 9, from which it is straightforward 
to determine the T*^*) that have energy-eigenvalues that depend upon a given phase-shift 5i ^. 
A cursory study of table [ijshows that ^^-states will, in general, receive contributions to their 
energy from interactions with / = 0, 4, 6, 8, as is well known [6j, and similarly for the other 
r*^*). As the dimensionality of an S0(3) irrep (which is 2Z + 1 for |/, m)) must be equal to the 
sum of the dimensionalities of the cubic irreps in its decomposition, cubic irreps will, in gen- 
eral, appear multiple times (with multiplicities denoted by N{T^^\l)) in the decomposition 
of an SO (3) irrep. Multiplicities greater than one occur for / > 5. The space associated with 
the j*^ occurrence of F^*-* in the decomposition of |/, m) is spanned by the orthonormal basis 
{ |rW^ Lz'-i }, where the number of values of equals the dimensionality of F*^*\ e.g. for 
/ = 5, the 3-dimensional irrep Tf occurs twice, and the space associated with the second 
occurrence is spanned by { |Tj~,0;5;2) , \T^ ,1;5;2) , |Tf,3;5;2)} ^. When calculating 
observables in a cubic volume, operators transforming as a component of a spherical tensor 
of rank-S, Og^\ are most conveniently written as 

= Yl 0^^'''^'^-^'^^ |F»,L,;5;j)(F«,L.;^;j| , (1) 
where the values of the are simply determined by matrix elements of Og 



^ The irreps of the full cubic group are F'^'^ = A^, yl^, E^^T^ , and T^, and have dimensionality 1, 1, 2, 3 
and 3 respectively. The superscript denotes the parity of F^^^. 

3 Each F(*) appears at least once in the decomposition of the \l,rn) with / < 6 except A-^ which first 
appears in the decomposition of the I — 9 irrep [SJj . It is important to note that the decompositions of 
the / = 7 and I — 8 irreps contain only F'-*-' that also appear in the decomposition of the / < 6 irreps, and 
consequently there is no F^'^ for which the dominant contribution to the interaction energy (in the large 
volume limit) is from the I — 7 and I = 8 partial- waves. 

^ The Lz quantum number indicates that a phase of e'^^''^ results from a (cubic) rotation oi cj) = mr/2 about 
the z-axis, with n an integer. = 3 is equivalent to a = — 1 and = 2 is equivalent to Lz — —2. 
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between \l,m), or |rW^ L^; S; j), or any states forming a basis in which the projections onto 
|r*^*\ Lz'-i S;j) are known. In determining the energy-eigenvalues of the states in the volume, 
it is the scattering amplitude in a given partial- wave that is written in the form of eq. ([T]), 
with S = I. 

The relations between the energy-eigenvalues of two-hadrons in a cubic volume and their 
scattering phase-shifts below the inelastic threshold, originally derived in the context of 
non-relativistic quantum mechanics, were shown to be valid in quantum field theory (QFT) 
without modification by Liischer [H The energy-shifts of scattering states due to the 
interactions exhibit power-law dependence upon the volume when the range of the inter- 
action is negligible compared to the spatial-extent of the volume. Corrections arising from 
the range of the interaction (for the case of Tr+Tr"*" the range is set by i? ~ l/(2m7r), while 
for nucleon-nucleon interactions it is set by i? ~ I/'^tt) are exponentially suppressed for 
L ^ R, and of the form ~ ^-^/R [55]. in this work, it is assumed that these finite-range cor- 
rections are negligible compared to the power-law energy-shifts due to the interactions. It is 
straightforward to calculate a two-hadron Green function resulting from an arbitrary source 
and sink. The Green-function is generated by the bubble-diagrams with non-interacting 
two-hadron states propagating from the source through multiple insertions of the T-matrix, 
and then to the sink. In free-space, the Green function exhibits poles at the location of 
bound-states and cuts along the positive real axis. In the finite- volume, modifications to 
the propagation of the two non-interacting hadrons eliminates the cuts on the positive real 
axis, replacing them with poles at the location of the energy-eigenstates. Further, these 
modifications shift the location of the poles on the negative real axis (if present in infinite 
volume). The energy-eigenvalues, corresponding to both bound-states and continuum states 
in the infinite- volume limit are determined by solutions to [8J, 

det cos 5 — sin 5 ] = , (2) 

where cot 5, sin 5 and F^^^'' are (/™^^ + 1)^ x (Z™*^^ + 1)2 dimensional matrices when the 
phase-shifts 61 are non-zero for / < Z™^^ and vanish for / > /™^^. Initially, it is convenient 
to work in the \l,m) basis in which, for uncoupled channels, cot 5 and sin 5 are diagonal 
matrices of the form 

cos 5 = cos 61^^12 Smum2 , siu 5 = sin 5i^^i^ ^mi.ma , (3) 

for h,2 ^ l™^^., but in which F^^^^ has off-diagonal elements, in general. F^^'^^ is a matrix 
that is a function of the dimensionless quantity ^ = f^, where q is related to the energy of 



the interacting two-hadron state, Eh^h2 = y Q'^ + "^l^i + y Q'^ + "^Hj' matrix elements 
are of the form 



^^Lm. = V(2/i + 1)(2/2 + 1) 



1712 

q TT^ 



l=\h-l2\ 

^ 1 { hi l2 \ { h 1 h 



E 



_ \000/ \ —mi — m m2 

m=—l 



ZjM-rq') , (4) 



where the functions Zi^^i)-] (f) are those defined by Liischer [H [2], 

2i,m{s;q ) = 2^ r \>2_z2 r , (5) 
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where Yim{^) are the spherical harmonics. The function ^^o.oll;?^) is UV-divergent and is 
defined with the same renormahzation scheme used to define the infinite-volume scattering 

amplitude. It is useful to first diagonalize the blocks of F'^^^^ with h = h = ^ F[^i\m2^ which 
diagonalizes F^^^^ into blocks with dimensions dictated by the number of occurrences of each 
r*^*) for / < Z™*", while leaving the sin 5 and cos 5 matrices diagonal. Further diagonalizations 
that may be required are confined within each r^*-*. The determinant in eq. (|2| becomes the 
product of determinants resulting from each V^^\ and therefore, for Z™*^^ = 6, this procedure 
requires dealing with matrices of size 4 x 4 or smaller. 

Despite the fact that eq. ([2]) requires forming the determinant of a finite dimensional ma- 
trix, it has infinitely many solutions. It is derived from a Green-function between arbitrary 
sources and sinks which, in principle, can couple to all of the eigenstates in the volume, 
manifested in the infinite-sums over integer-triplets that define the ^;,m-functions. There- 
fore, the zero's of the determinant in eq. g define all of the energy-eigenvalues and hence 
eigenstates. As discussed previously, the energy-spectrum of two non-interacting hadrons in 
the cubic volume with periodic BC's, and with vanishing total momentum can be defined 
by triplets of integers, n. 




E = Y |qip + ml + Y |q2p + ml \l { — ] l^p + rnf + \l { — ] |np + ml 

where one hadron carries momentum qi = and the other carries momentum q2 = — ^n, 
and the reduced mass of the system is = m^^ + m^^. This (non-interacting) spectrum 
is recovered in the above formalism, in particular eq. ([2]), in the limit that 5/ — )■ in each 
partial-wave from the poles in the ^;_m-f unctions that exist along the positive real axis. 
The degeneracy of any given |np-shell is straightforward to determine and is recovered 
from the number of states in the r*^*-* that span the |np-shell, as shown in table |ll| As the 
(single hadron) momentum eigenstates in a given |np-shell are degenerate, the corresponding 
r*^*) are also degenerate. These degeneracies are lifted by two-particle interactions that 
induce non-zero 6i^s. Table [n] shows that all but one of the F^*) are required to describe the 
eigenstates for |np < 6, and from table [l] it can be concluded that for 5; 7^ for / < 6 all 
of the eigenstates with |np < 6 are shifted from the non-interacting two-hadron energy due 
to interactions. However, the A~[ irrep first occurs in the |np = 14 shell and its energy is 
dependent upon interactions with / > 9. 



III. ENERGY-EIGENVALUES, SOURCES AND SINKS 

Liischers formalism, as detailed in the previous section, is used to construct explicit relations 
between the energy-eigenvalues of the F*^*^ and the interaction phase-shifts for I < 1^°"^ = 6, 
the results of which are presented in this section. Sources and sinks for LQCD calculations 
that transform as a given F^*-* are constructed from the single-hadron momentum-eigenstates, 
and Fourier transformed into position- space. One pair of these sources and sinks would 
couple only to a single energy-eigenstate in the absence of interactions between the hadrons. 
As the interactions do not induce mixing between distinct T^'^\ these sources and sinks 
couple, in principle, to all states that transform in the same F*^*-*. To keep the presentation 
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TABLE II: The degeneracies of, and the number of occurrences of each T^*) in, the lowest-lying 
|np-shells. Note: the irrep first appears in the |np = 14 shell. 
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of results simple, explicit derivations are deferred to Appendix [A], where calculations of 
the even- and odd-parity systems with Z™'"^ = 4 are detailed, and which straightforwardly 
generalize to any 

As the hadronic interactions considered in this work result from QCD with the strong CP- 
violating parameter 6 set equal to zero, and without the electroweak interactions, parity is a 
good quantum number. Consequently, the contributions to the finite- volume function F^^'^^ 
do not mix states of opposite parity, and therefore the required calculations decompose into 
the parity-even and parity-odd sectors. If weak-interactions are included in the analysis, as 
will necessarily be the case when hadronic parity-violating interactions are calculated with 
LQCD, mixing between the parity-sectors will occur. 



A. Positive Parity Systems 



There are five positive parity irreps of the cubic group, A^, A2, E^, T^, and with 
dimensions 1, 1, 2, 3, and 3 respectively. Table|T]shows how the interactions in a given partial- 
wave contribute to each F^*^. The energy-eigenvalues, sources and sinks for the even-parity 
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IIIA2, 
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III A 4 


and 


III A 5 



1. A'l Representation 

The energy-eigenvalues of A~i^ states depend upon the phase-shifts in the / = 0, 4, 6, 8, ... 
partial-waves, as can be seen in table |T] Diagonalization of the blocks in the finite-volume 
function of the form F^f^^ for / = 0, 4, 6 gives the states \ A'f, Lz] I; j), as defined immediately 
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before eq. ([T]), with 

|A+,0;0;1) 

0;4;1) 
|A+ 0;6;1) 



|0,0) 
4 



6,4) 



2^2 



|6,0) + 



4 



|4,-4) 
4) 



(7) 



for the Al eigenstate of each ' ' in the orbital angular momentum (spherical-wave) basis 

|Z,m). With these states and the corresponding eigenvalues from f\^j^\ the procedures 
described in Appendix |A] allow for the contribution to eq. ^ from states to be written 
as 

/p(Fy,A+) 

-T n 



det 



^cot(5o 

COt54 

cot^fi 



0;0 

{FVAt) 
4;0 

\^6;0 



0;4 
(FVAt) 
4;4 

(FVAt) 
6;4 



0;6 
4;6 

-.{FVAt) 







(8) 



where the finite-volume contributions are 
-^{FVAt) ^0,0 (l;*?^) 
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0:4 



-^{FVAt) 
^ 0:6 
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1008 
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;;2> 



(FVAt) _ 
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and fIJ^^'^') 



2016v/2lZio,o(l;g^) 30492Zi2,o (1; 
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371457r3/2gi3 
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. Equation (8) yields an infinite number of energy-eigenvalues and 
eigenstates, each of which depend upon the phase-shift in the / =0, 4, and 6 partial waves. 

In the |np-shells for which there is just one A\ state, as shown in table|ll| its energy-shift 
due to interactions receives contributions from the I = 0,4,6, ... partial-waves. However, in 
the |np-shells in which there are multiple states (first occurring at |np = 9), the energy- 
eigenstates are linear combinations of these states. In the large- volume limit, the shift in 
the energy-eigenvalue of one combination is dominated by the interactions in the I = 
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FIG. 1: The function -Zq.o (ij 9^)- The vertical dashed lines denote the position of the poles of the 
function corresponding to the non-interacting energy-eigenvalues. 



partial-wave, the shift in a second combination is dominated by the interactions in the I = 4 
partial-wave, the shift in a third combination is dominated by the interactions in the / = 6 
partial-wave, and so on. So while the naive argument that Af states receive contributions 
from interactions in the / = 0,4,6,... partial-waves is generally true, linear combinations 
of Af states are formed such that it is not true in the infinite- volume limit. The energy- 
shift of each occurrence of an A^ energy-eigenstate in a given |np-shell is dominated by 
the interaction in a different partial- wave in the infinite- volume limit. To demonstrate this 
point, consider the situation where the phase-shift in the / = 6 partial- wave vanishes, in 
which case eq. (|8| becomes a 2 x 2 matrix with the following two solutions: 

cot^o , cot(54 ^0,0 (!;<?') 280^8,0 40^6,0 54^4,0 (I;?') 



3/2, 



q 143yT77r3/2g9 Ilv^7r3/2g7 USiT^/^q^ 



2 V V143v^7r3/V llv/T37r3/V U3n^/2q5 



= 0. (9) 



In the case of tan54 ^ tanSo, the / = dominated solution is 
qcotSo 



2 .2 122887r^Z4,o(l;g^)^ . , . 

Zo,o{l;q) + ^YTo r + O (tan 54) , (10) 



^/KL 



7L10 [g9cot54 
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and is valid for all |np-shells. If phase shifts in both the / = 4 and / = 6 partial- waves 
vanish, eq. rt8| and eq. (10) reduce to the familiar result found by Liischer ^, 



gcotf^o 



(11) 



where the function Zq^ (1; g^) is shown in fig. [Tj Performing a large-volume expansion of 
the solution (as discussed in Appendix [C|) to eq. (11) in the |np = 9-shell gives the energy- 
eigenvalue 



1 

2/1 



367r2 20tan5o(|n|2 = 9) 



+ ^(tan^^o) + 



L2 L2 

while the second solution to eq. ^ has a perturbative expansion of the form 



(12) 



1 

2^ 



367r2 8960tan54(|n|2 = 9) 



L2 



243 L2 



+ ^(tan^e) + ••• 



(13) 



where the contribution from the / = partial-wave is strongly suppressed in the large- 
volume limit. While the two basis states, ||n|2 = 9; Af{l)) and ||n|2 = 9; Af{2)), both have 
a non-vanishing overlap with \l,m) = |0,0), it is obvious that a linear combination can be 
formed that has vanishing overlap. Inserting the interactions once, as is appropriate for 
determining the energy-eigenvalues in large volumes {i.e. first order perturbation theory in 



1/L), dictates the form of the expansions in eq. (12) and eq. (13). 

Sources and sinks that have an overlap, and in general a range of overlaps, with the 
finite-volume energy-eigenstates of hadronic systems are required for LQCD calculations. 
While the interactions between hadrons gives rise to energy-eigenstates that are not products 
of single-hadron eigenstates of the linear-momentum operator, sources and sinks can be 
constructed from the single-hadron momentum eigenstates that transform as a given r*^*\ 
that will have non-zero overlap with the energy-eigenstates in the same irrep. Constructing 
sources and sinks from single hadrons that have equal and opposite momenta ensures that 
the total momentum of the combined system vanishes. The relative-momentum-eigenstates 
of definite parity, V, are denoted by 



|n , V) 



where V is the parity of the state {V = ±1) and n 




(n^O) 

(n = and V 



+ 1) 



(14) 



is the triplet of integers that 



define the relative momentum of the two-body system. The states in eq. ( 14 ) are eigenstates 
of the relative kinetic energy operator Trei, with the eigenvalues displayed in eq. ([6]). By 
taking appropriate linear combinations of these momentum-eigenstates, states in the 
representation (or any other irrep) can be constructed in each |n|2-shell if the shell supports 
it (see table [n] and Ref. [53]). For example, in the |np = shell the basis-state is 



|n|2 = 0; At) = 1(0,0,0) , V = +1) 



^ The "S-function" , S{q'^), used in, for example, Ref. [S5], is related to -Zo,o(l; 9^) by S'(g^) = \/47rZo,o(l; q^)- 
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TABLE III: The momentum-space structure of Af sources and sinks for |np=0-3. These are 
shown graphicahy in figs. [2] and |3j 



n 



1(0,0,0) , +1) 1 



n 



1(1,0,0) , +1) 
1(0,1,0) , +1) 
1(0,0,1) , +1) 



1 
1 

V3 

1 

V3 



|np=2 
1(1,1,0) , +1) 
1(1,0,1) , +1) 
1(1,0,-1) , +1) 
1(1,-1,0) , +1) 
1(0,1,1) , +1) 
1(0,1,-1) , +1) 



1 

V6 
1 

1 

V6 
1 

Ve 
1 

V6 
1 

V6 



n 



1(1,1,1), +1) 
1(1,1,-1) , +1) 
1(1,-1,1) , +1) 
1(1,-1,-1) ,+1) 



TABLE IV: The momentum-space structure of Af sources and sinks for |np=4-6. These are 
shown graphicahy in figs. [3] and |4j 



n 



n 



n 



1(2,0,0) ,+1) ^ 
1(0,2,0) , +1) ^ 
1(0,0,2) ,+1) ^ 



1(2,1,0) , +1) 
1(2,0,-1) , +1) 
1(1,2,0) , +1) 



2^3 
1 

2^3 



1(1,0,-2) , +1^ 1 



2^3 
1 

1(0,2,1), +1) ^ 



1(2,0,1), +1) ^ 

1(2,-1,0) , +1) ^ 

1(1,0,2), +1) ^ 

1(1,-2,0) , +1) ^ 

1(0,2,-1) , +1) ^ 



1(2,1,1), +1) ^ 

(2,-1,1) , +1) ^ 

1(1,2,1) , +1) ^ 

1(1,1,2), +1) ^ 

(1,-1,2) , +1) ^ 



1(2,1,-1) , +1) ^ 

l(2rl,-l) , +1) ^ 

1(1,2,-1) , +1) ^ 

1(1,1,-2) , +1) ^ 

1(1,-1,-2), +1) ^ 

l(W),+i) ^ 



while for |np = 1, the basis-state is 

1(1,0,0) ,V = +1) + 1(0,1,0) ,V = +1) + 1(0,0,1) , P = +l) 



n 



1; At) 



In general the coefficients of these basis vectors are valid up to an arbitrary phase. The 
complete momentum-space basis for the Af sources and sinks in each |np-shell are presented 
in table III for |np < 6. 

The momentum-space representations of the sources and sinks (the left panels in figs. [2]- 
4]) show the n- vectors that transform as an Af in the given |np-shells. The widths of the 
vectors are proportional to the magnitude of their amplitudes and their color denotes the 
sign (red=positive, blue=negative). The position-space representations of the sources and 
sinks (the right panels in figs. |2||4|) show the surfaces of constant pn,p(i'), defined by 



Pn,p(r) = K r I {n^,ny,n^) , P )P , (15) 

which are obtained by Fourier transform. In the position-space representations, r refers to 
the relative distance between the two particles and in the figures r is defined to be r = r/L. 



12 



A[ : |Hp =0 : = 




At : |rtp = 1 : L- = 




: |Hp = 2 : L- = 




FIG. 2: The momentum-space representations (left) and position-space representations (right) of 
two-body relative states in the Af representation for select |np shells. Here f = r/L. 
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A[ : |Hp = 3 : = 




At : |rtp =4 : L- = 




: |«p = 5 : L- = 




FIG. 3: The momentum-space representations (left) and position-space representations (right) of 
two-body relative states in the Af representation for select |np-shells. 
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FIG. 4: The momentum-space representations (left) and position-space representations (right) of 
two-body relative states in the representation for |np=6 shell. 



2. A2 Representation 

The other one-dimensional positive-parity irrep of the cubic group is the A2 . Due to its 
complexity, the lowest-lying state transforming as a A2 is in the |np = 5 shell, as indicated 
in table [TTj Further, the lowest partial-wave contributing to its energy is / = 6, and as this 
analysis is truncated to partial- waves with / < 6, the contribution to the determinant in 
eq. ([2]) has the solution 



g^^cot^e 



L J 7r3/2 

6g8Z4,o (1; e) 160^g6^6,o (1; q') m^Z^fi (1; 



17 323 19^ 

2592x/2Tg2Zio,o (1; q") ^ 1980Zi2,o (1; g^) ^ 264v/T00TZi2,4 (1; . 



2nV' 1 



7429 7429 7429 

m ■ (16) 



L J 7r3/2 

and the associated eigenstate of the F^^^^^ block is 

\Ai,2:6;l} = \^\e,2) + ^^116.-2) - - jyf|6.-6> . (17) 

The function is shown in fig. [s] as a function of g^. Its pole at g^=5, denoted by the 
vertical dashed line, corresponds to the non-interacting {Sq = 0) energy-eigenvalue. This 
is the only |np-shell with |np < 6 which supports the A2 irrep, as shown in table [nj In 
fig. [6] we give the graphical representations of the source and sink that generates this irrep 
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2e+06 



1 .5e+06 



1e+06 



500000 



-500000 



-1e+06 



-1.5e+06 



-2e+06 




2 

~2 



FIG. 5: The function Xj^^ , as defined in eq. ( 16 ) , as a function of . The vertical dashed hne denotes 
the position of the pole of the function corresponding to the non-interacting energy-eigenvalues. 



TABLE V: The momentum-space structure of the source and sink in the |np 
are shown graphically in fig. [6} 



5-shell. They 




1(2,0,1) , +1) 
1(2,-1,0) , +1) 
1(1,0,2) , +1) 
1(1,-2,0) , +1) 
1(0,2,-1) , +1) 
1(0,1,-2) , +1) 



1 

]_ 

]_ 

2V3 
1 

2^/3 

]_ 

2V3 
1 

2V3 



in the |np = 5-shell. As this is the lowest-lying state whose energy-eigenvalue is insensitive 
to / < 6 interactions, it is LQCD correlation functions constructed to transform in the 
irrep that will enable a calculation of 6q. However, as the lowest energy contributing to an 
A2 correlation function occurs in the jn^l = 5-shell, relatively large lattice volumes will be 
required in order to have this state lie below the inelastic threshold. The momentum-space 
structure of the source and sink that couple to the A2 state in the |np = 5-shell is given in 
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FIG. 6: The momentum-space representation (left) and position-space representation (right) of 
two-body relative states in the representation for the |np = 5 shell. 

table lYl 

3. Representation 



The energy-eigenvalues of states transforming in the irrep receive contributions from 
interactions in the / = 2,4,6,... partial-waves. As the irrep is two-dimensional, the 
contribution to the determinant in eq. ^ results from a 6 x 6 matrix when / < 6. However, as 
the two states in the E^ irrep (with = and Lz = 2) are degenerate, and orbital-angular 
momentum is conserved by the interactions (unlike the situation in the baryon-sector) , the 
analysis can be reduced to that of a 3 x 3 matrix. The E^ Lz = states associated with 



the F. 



(FV) ^(FV) 
2;2 



4:4 



, and F^^Q^ blocks are 



|E+,0;2;1) 
|E+,0;4;1) 

|E+,0;6;1) 
and the contribution to eq. ([2]) becomes 



|2,0) 
1 



4,4) - y^|4,0) + y!|4,-4) 



(18) 



4,6,4) + y^|6,0) + i|6,-4) 



det 



'cot52 
cot(54 
cot^g/ 



/^{FV,E+) -^{FV,E+) -^{FV,E+)\ 
2-2 2-4 2 6 

iFV,E+) -Y^{FV,E+) -y^{'fV,E+) 



-^4;2 -^4;4 

{FV,E+) -r={FV,E+) Y^{FV,E+) 



-^4;P 



6:2 



6:4 



6;6 



(19) 
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where 

-=;(Fy,E+) 



2;2 



^{FV,E+) 
' 2;6 



-^{FV,E+) 



F, 



{FV,E+) 
4;6 

iFV,E+) 
6;6 



7]-3/2g 7^3/2g5 



30^^^6,0(1; g') 40v^Z4,o (1; g') 



Il7r3/2g7 



+ 



777r3/2g5 

4./fZ6,o(l;g 



30. 



^^4,0 (1; 9^) 



Il7r3/2g 



7r3/2g9 Il7r3/2g7 ' Il7r3/2g5 

Zo,o(l;g^) 392^8,0 (l;g^) 64^6,0 (l;g^) 108^4,0 (1; g^) 
^3/2;; 143v^7r3/2g9 Ilv^7r3/2g7 100l7r3/2g5 

128^^0 (1; e) 18v/210Z6,o (1; g') (1; Q 

Il7r3/2g5 



7r3/2g 

1512^/^2^10,0 



3237r3/2gll 2097r3/2g9 

Zo,o(l;g^) ^ 30492Zi2,o(l;9 , ^ 



7r3/2g 



371457r3/2g 



ol3 



371457r3/2g 



;513 



1877r3/2g7 

+ 



. ;^2-i 



264^l00lZi2,4(l;g') , 1152y2lZio,o (1; 



+ 



280^8,0 (l;g2) 480v/l3Z6,o (1; g^) 114^4,0 (1; 

209VT77r3/2g9 35537r3/2g7 1877r3/2g5 



74297r3/2gii 
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TABLE VI: The momentum-space structure of E^, Lz = sources and sinks for |n|2=l,2, and 4. 



n 



1(1,0,0) , +1) 
1(0,1,0) , +1) 

1(0,0,1) , +1) 



1 

Ve 
1 



n 



1(1,1,0) , +1) 
1(1,0,1) , +1) 
1(1,0,-1) , +1) 
1(1,-1,0) , +1) 
1(0,1,1) , +1) 
1(0,1,-1) , +1) 



1 

V3 

1 
1 

2^3 
1 

V3 
1 

2^3 
1 

2v^ 



n 



1(2,0,0) , +1) ^ 
1(0,2,0) , +1) ^ 
1(0,0,2) , +1) - 



It is obvious that the solutions of eq. (19) depend upon the / = 2, 4, and 6 partial-waves in 
a non-trivial manner. 

In the limit of vanishing interactions in partial-waves with / > 4, the contribution from 

;su 

Ao(i;g') 



the irrep to eq. ( 19 ) results from a 2 x 2 matrix, and has solutions 

69^4,0 (l;g') 



COt52 _|_ COt54 



7]-3/2g 



z2\ 



196^8,0 (l;g2) 32^6,0 (1; 

143v^7r3/2g9 Ilv^7r3/2g7 1437r3/2g5 



±- 



1 



/392Z8,o(l;g') 64^6,0 (l;g') 750^4,0 (1; g'^ 



Vl43yT77r3/2g9 llv/T37r3/2g7 100l7r3/2g5 

30 



+ COt(52 — COt54 



+4 



21 1/2 



13^6,0 (1; g ) ^ 40^3^4,0 (1; g^) 



Ilvr3/2g7 



777r3/2g5 



(20) 



In the limit that tan54 << tan52, the / = 2 dominated solutions to eq. (|20|) result from 
g^cot52 = 



27r 



L 



TT 



3/2 



6 



g%o(i;g') + ^Ao (i;g') 



27r 



TT 



3/2 



^1 (g^) , (21) 



where function is shown in fig. [Tjas a function g^ ^. The graphical representations of 
the sources and sinks that generate this irrep for |n|2 < 6-shell are shown in figs, [s] and [9] in 
the case of Lz = 0, and the momentum-space structures are given explicitly in table |VT] and 
table IVnH 

There are two occurrences of the irrep in the jnp = 5-shell. Linear combinations of 
the basis states can be formed: one that is dominated by ^2, and one that is dominated by 



6a in the infinite- volume limit, as shown in table VI and table VIII The states are defined 



by {l,m\\n\f;T^'^\Lz) = (2, 0|54; -E'^, 0) = 0, and the orthogonal combination |52;-E^,0). 
As is the case in the sector, these states are not energy-eigenstates since they have a 
non-zero projection, in principle, onto all states. The perturbative expansions of the 
energy-eigenvalues in the large- volume limit can be found in Appendix [C| 



This expression has been derived previously by R. Briceno |19) . 
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: |Hp = 1 : = 2 




: |rtp = 2 : L- = 2 




: |«p =4 : L- = 2 




FIG. 8: The momentum-space representations (left) and position-space representations (right) of 
two-body relative states in the representation with = 2 in the |np=l, 2, and 4 shells. 
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£+ (1=2) : \nf = 5 : = 2 




£+(1=4) : \nf = 5 : L - = 2 




: |Hp = 6 : L- = 2 




FIG. 9: The momentum-space representations (left) and position-space representations (right) of 
two-body relative states in the representation with = 2 in the |np=5 and 6 shells. 
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TABLE VII: The momentum-space structure of E~^, = sources and sinks for |np=5 and 6. 



n 



\l=2) 



1(2,1,0) , 
1(2,0,1) , 
1(2,0,-1) , 
1(2,-1,0) , 
1(1,2,0) , 
1(1,0,2) , 
1(1,0,-2) , 
1(1,-2,0) , 
1(0,2,1) , 
1(0,2,-1) , 
1(0,1,2) , 
1(0,1,-2) , 



+1) 
+1) 
+1) 
+1) 
+1) 
+1) 
+1) 
+1) 
+1) 
+1) 
+1) 
+1) 



2V78 
1 

V78 
1 

5 



2V78 
5 



2V78 
7 



2V78 
5 



2V78 
1 

VT8 
1 

Vts 

7 

2V78 
7 



2V78 



|np=5(;=4) 



1(2,1,0) , 

1(2,0,1) , 
1(2,0,-1) , 
1(2,-1,0) , 
1(1,2,0) , 
1(1,0,2) , 
1(1,0,-2) , 
1(1,-2,0) , 
1(0,2,1) , 

1(0,2,-1) , 
1(0,1,2) , 
1(0,1,-2) , 



+1) 
+1) 
+1) 
+1) 
+1) 
+1) 
+1) 
+1) 
+1) 

+1) 
+1) 
+1) 



2^26, 
2_ 
13 

2_ 
13 



2V26 
3 



2V26 
1 



2V26 
1 



2V26 
3 



2V26_ 
2_ 
13 

2_ 
13 

1 



2V26 
1 



n 



=6 



1(2,1,1) 
1(2,1,-1) 
1(2,-1,1) 
1(2,-1,-1) 
1(1,2,1) 
1(1,2,-1) 
1(1,1,2) 
1(1,1,-2) 
1(1,-1,2) 
1(1,-1,-2) 
1(1,-2,1) 
1(1,-2,-1) 



,+1) 
,+l> 
,+l> 
,+1) 
,+1) 
,+1) 
,+1) 
,+1) 
,+1) 
,+1) 
,+1) 
,+1) 



1 

'2^6 
1 

'2^6 
1 

'205 
1 

2^6 
1 

'2^6 
1 

' 2^6 
1 

1 

V6 
1 

Ve 
1 

Ve 
1 

~ 2^/6 
1 

~ 2Ve 



TABLE VIII: The momentum-space structure of E^, = 2 sources and sinks for |np=l, 2, and 
4. These are shown graphically in fig. [8} 



|n|2=l 


|n|2=2 


|n|2=4 


1(1,0,0) , +1) 
1(0,1,0) , +1) ^ 


1(1,0,1) , +1) -i 
1(1,0,-1) , +1) 
1(0,1,1) , +1) 1 
1(0,1,-1) , +1) 1 


1(2,0,0) , +1) 
1(0,2,0) , +1) ^ 



4- Representation 

The energy-eigenvalues of states transforming in the irrep receive contributions from 
interactions in the / = 4, 6, ... partial-waves. The irrep is three-dimensional, with states 
identified by = 0, 1, 3, and provides a contribution to the determinant in eq. (g that 
results from a 6 x 6 matrix for / < 6. As the three L^-states are degenerate, the analysis 

collapses down to that of a 2 x 2 matrix. The = states associated with the F^^^^ 

and F^^^^ blocks are 

|T^i',0;4;l) = i=|4,4) - -^\^,-^) 

Y 1 (22) 

|T+,0;6;1) = -^|6,4) - ^|6,-4) , 
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TABLE IX: The momentum-space structure of 
These are shown graphically in fig. [9| 



2 sources and sinks for |np=5 and 6. 



n 



^(^=2) 



1(2,1,0) 
1(2,0,1) 

1(2,0,-1) 
1(2,-1,0) 
1(1,2,0) 
1(1,0,2) 
1(1,0,-2) 
1(1,-2,0) 

1(0,2,1) 
1(0,2,-1) 
1(0,1,2) 
1(0,1,-2) 



+1) 
+1) 

+1) 
+1) 
+1) 
+1) 
+1) 
+1) 

+1) 

+1) 
+1) 
+1) 



2V26 
2_ 
13 

2_ 
13 



2^26 
3 

2v^ 
1 



2V26 
1 



2^26 
3 



2V26 
_2_ 
13 

2_ 
13 



2V26 
1 



2V26 



n 



1(2,1,0) 
1(2,0,1) 
1(2,0,-1) 
1(2,-1,0) 
1(1,2,0) 
1(1,0,2) 
1(1,0,-2) 
1(1,-2,0) 
1(0,2,1) 
1(0,2,-1) 
1(0,1,2) 
1(0,1,-2) 



+1) 
+1) 
+1) 
+1) 
+1) 
+1) 
+1) 
+1) 
+1) 
+1) 
+1) 
+1) 



2V78 
1 

1 

5 



2V78 



2V78 
7 



2V78 
7 



2V78 
5 

2^78 

]_ 

V78 

]_ 

•/T8 
7 



2V78 
7 

2V78 



|nP=6 



1(2,1,1) , 
1(2,1,-1) 
1(2,-1,1) 
1(2,-1,-1) 
1(1,2,1) , 
1(1,2,-1) 
1(1,-2,1) 
1(1,-2,-1) 



+1) 
+1) 
+1) 
-+1) 
+1) 
+1) 
+1) 
>+l) 



]_ 

2^2 
]_ 

2V2 
]_ 

2V2 
]_ 

2V2 

1 

2V2 
1 

2V2 
1 

2V2 
1 

2V2 



and the contribution to eq. ^ is 



where 



4;4 



F, 



iFV,T+) 
4;6 

(FV,T+) 
6;6 



det 



COt54 

cot^e 



^(FV,T+) -^{FV,T^r 

4;4 -T 4-6 

^(FV,T+) p(Fy,T+) 

'-^6;2 -^6;6 



Zo,o(l;g') 448^8,0 (l;g') 4^6,0 (l;g=^) , 54^4,0 (l;g 



;^2\ 



7r3/2q, 



576 1 



'21 7 



3237r3/2gii 



143VT77r3/2g9 

(l;g2) 112 
- + 



llVT37r3/2g7 

Ao(i;g') 

- + 



1437r3/2g5 

42^5^6,0 (1; e) 



2097r3/2g' 



1877r3/2g 



;77 



(23) 



12^^^4,0(1; g' 

Il7r3/2g5 



;;2\ 



Zo,o (1; g') 26136Zi2,o (1; <?') ^ 1584v/T00TZi2,4 (1; g') , 624v^Zio,o (1; g') 



7]-3/2g 



+ 



371457r3/2gl3 ' 371457r3/2gl3 

120^8,0 (l;g^) 80v^Z6,o (1; g^) 96^4,0 (l;g^) 
209v^7r3/2g9 35537r3/2g7 1877r3/2g5 



74297r3/2gll 



23 



80000 



60000 



40000 



20000 



-20000 



-40000 



-60000 



-80000 
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FIG. 10: The function Af^, as defined in eq. (25), as a function of (p. The vertical dashed 
hnes denote the position of the poles of the function corresponding to the non-interacting energy- 
eigenvalues. 



The solutions to eq. (23) are obtained from 



cotSi ^ cot^e 312721^10,0 (1; q^) ^ 13068Zi2,o (1; q^) 



±- 



+ ■ 



792v/TooTZi2,4(l;g^^ 



74297r3/2gii 



371457r3/2gl3 



Zo,o(l;a ^ 15^4,0 ^ 106^6,0 (1; g^) ^ 



371457r3/2gi3 



624V21Z 



'10,0 



22l7r3/2g5 ■ 323v^7r3/2g7 ' 247v^7r3/2g9 

26136Zi2,o(l;g') 



1584v/l00lZi2,4(l;g^) 2166^4,0 



74297r3/2gll 371457r3/2gl3 

252^6,0 



371457r3/2gl3 



243l7r3/2g5 



576 1 



+ 



21 7 

65^10,0 



^ 10072^8,0 



3553v^7r3/2g7 2717v^7r3/2g9 

12./^Z4n(l;g^) 



COt54 — COt5f 



323g 



nil 



42v^Z6,o(l;g'^ 



112 



221^8,0 (1; Q'^^ 



187g 



209g 



24 




FIG. 11: The momentum-space representations (left) and position-space representations (right) of 
two-body relative states in the representation with = in the |np=5 and 6 shells. 



In the situation where the interaction in the / = 6 (and higher) partial-wave vanishes, the 
energy-eigenvalues are sensitive to the I = 4 interaction alone, and can be found from 



1 

TJ ^ 

Zn 

T 



g%o(i;g') 



448^8,0 (l;g') 4g2Z6,o(l;g 



143V17 



11V13 



143 



3/2 



VT' 



(25) 



10 



where the function (g^) is shown in fig 

The graphical representations of the sources and sinks that generate the irrep for the 
low- lying |np-shells are shown in fig. 11 
given explicitly in table |X], and for = 
and sinks are related to those with L, - 



and the momentum-space structures for = are 



1 in table XI The structures of the = 3 sources 



1 by complex conjugation of the coefficients. 
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TABLE X: The momentum-space structure of , = sources and sinks. These are shown 



graphicahy in fig. 11 



|np=5 


|np=6 


1(2,1,0) , +1) 1 
1(2,-1,0) , +1) -\ 
1(1,2,0) , +1) -i 
1(1,-2,0) , +1) \ 


1(2,-1,1) , +1) 
1(1,2,1),+!) 
1(1,-2,1) , +1) ^ 


1(2,1,-1) , +1) ^ 
1(2,-1,-1) , +1) 
1(1,2,-1) , +1) 
1(1,-2,-1) , +1) ^ 



TABLE XI: The momentum-space structure of , Lz = 1 sources and sinks. 



n 



n ■ 



=6 



1(2,0,1) 
1(1,0,2) 
1(0,2,1) 
1(0,1,2) 



+1) 
+1) 
+1) 
+1) 



1 

2^/2 
]_ 

i 

2V2 
_ i 
2V2 



1(2,0,-1) 
1(1,0,-2) 
1(0,2,-1) 
1(0,1,-2) 



+1) 
+1) 
+1) 
+1) 



i_ 

2V2 
1 

2V2 

_ i 
2V2 
i 

2V2 



1(2,1,1) 
1(2,-1,1) 
1(1,2,1) 
1(1,1,2) 
1(1,-1,2) 
1(1,-2,1) 



+1) 
+1) 
+1) 
+1) 
+1) 
+1) 



- + - 



1(2,1,-1) 
1(2,-1,-1) 
1(1,2,-1) 
1(1,1,-2) 
1(1,-1,-2) 
1(1,-2,-1) 



+1) 
+1) 
+1) 
+1) 
+1) 
+1) 



- + - 
4^4 

1 _ i 

4 4 



5. T2 Representation 

The energy-eigenvalues of states transforming in the irrep receive contributions from 
interactions in the I = 2,4,6,... partial-waves. The irrep is three-dimensional, with 
states defined by = 1, 2, 3, and provides a contribution to the determinant in eq. Q that 
results from a 12 x 12 matrix for / < 6 (it is 12 x 12 and not 9x9 because there are two 
T2^s in the decomposition of / = 6, see table [l|. As the three L^-states are degenerate, the 
analysis collapses down to that of a 4 x 4 matrix. The = 2 states associated with the 

Fif/\ and Ff/^ blocks are 

17^2+, 2; 2; 1) = ^(|2,2) - |2, -2)) 



17^2*", 2; 4; 1) = ^(|4,2) - |4,-2)) 
17^2^2; 6; 1) = ^(|6,2) - |6, -2)) 
17^2^, 2; 6; 2) = ^(|6,6) - |6, -6)) 



(26) 
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FIG. 12: The function 'Vj^, as defined in eq. (30), as a function of (p. The vertical dashed fines 
denote the position of the poles of the function corresponding to the non-interacting energy- 
eigenvalues. 



With these basis states, the contribution to the determinant in eq. ([2]) becomes 

-,(FV,T+) 

^{FV,T+) 



det 



/ cot 5' 



v 



2 

cot54 

cot^e 

cot^e/ 



4;2 

{FV,T+) 
6i;2 
. {FV,T+) 
V 62;2 



{FV,T+) 
2;4 

{FV,T+) 
4;4 

{FV,T+) 
6i;4 

iFV,T+) 
62;4 



^ 2;6i 
-piFV,T+) 
^ 4;6i 

-^6i;6i 

-^62;6i 



-^(FV,T+)\ 
^ 2;62 

^ 4;62 

-^61,62 

^ 62,62 



(27) 
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where 

iFv,T+) Zofi (1; g') 4^4,0 (1; q^) 



F. 



2,2 ^3/2q 77r3/2g5 

{Fv,T+) 40yjz6,o(l;g^) 20^3^4,0 



F 



2,4 Il7r3/2g7 777r3/2g5 

(i.y,T+) _ Zo,o (1; q') 54^4,0 (1; g') , 20^6,0 (1; g') 



4,4 



7r3/2g 777r3/2g5 1 1 v/l37r3/2g7 

(FF,T+) _ 5-^^^4,0(1; g') 5v^Z6,o(l;g') 

2,6i - Ilvr3/2g5 Il7r3/2g7 

42-2'_n.^2\ in. . n. ;^2^ 



+ 28, /^Z8,o (1; q') 10 J^Ao (1; g') Vf ^6,0 (i; g') 1008, /^Zio,o (i; g' 

i^4R; ' = — ^ — + — — + V 



4,61 197r3/2g9 ^ Il7r3/2g5 ^ 1877r3/2g7 ^ 3237r3/2gll 



p(i^y,T+) _ 45^8,0 (l;g^) (1; g') 59^4,0 (1; g') , 620713^6,0 (1; g') 



19^7r3/2g9 7r3/2g 1877r3/2g5 ^ 35537r3/2g7 



^ 162v/2lZio,o(l;g^) ^ 3267Zi2,o (1; g') 19871001^12,4 (1; g') 



74297r3/2gll 74297r3/2gl3 74297r3/2gl3 



(Fv,T+) Ao (1; g^) y i^6,o (i; g^) 64^ j^^s.o (i; g^) 

-Tofi. = tttt;;^ h 



2,62 ^3/2^5 ^ 7,3/2^7 37^3/2^ 



^9 



^{FV,T+) 



yFv,T+) '^\J iMi^4,o (1; g^) Sy'^^Ze^o (1; g^) 20^ ^yZg^o (1; g^) 336^/^Zio,o (1; g^'' 



4,62 ^3/2^5 177r3/2g7 ^ 197r3/2g9 ^ 3237r3/2gll 



_^j,y^^+^ _ (1; g') I40^ff ^6,0 (1; g') 5^/ J^Zg^o (i; g') 666^/^Zio,o (i; g'^ 



61,62 177r3/2g5 ^ 3237r3/2g7 ^ 577r3/2g9 74297r3/2gll 

1287,/^Zi2,o (1; 858,/f Zi2,4 (1; g') 



74297r3/2gl3 74297r3/2gl3 



(FF,T+) _ Zo,o(l;g^) 9^4,0 (l;g^) _ 20Vl3Z6,o (1; g^) 5^8,0 (1; if) 18V2l2io,o (1; g') 



^^'^^ 7r3/2g ^ 177r3/2g5 3237r3/2g7 ^ 19v^7r3/2g9 74297r3/2gll 

23991Zi2,o(l;g') _ 594v/l00lZi2,4 (1; g') 

371457r3/2gl3 371457r3/2gl3 

(28) 



The solutions to eq. (27) must be determined numerically and will, in general, depend upon 
the interactions in the Z = 2, 4, and 6 partial-waves. In the limit where the interactions in 
the / = 6 and higher partial-waves vanish, leaving contributions only from interactions in 
the / = 2,4 partial- waves, the contribution to the determinant in eq. (27) collapses down to 
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TABLE XII: The momentum-space structure of T2 , = 2 sources and sinks for |np=2, 3, and 



5. These are shown graphicahy in figs. 13 



|n|2=2 


|n|2=3 


|n|2=5 


1(1,1,0) , +1) 
1(1,-1,0) , +1) ^ 


1(1,1,1),+!) \ 
1(1,1,-1) , +1) \ 
1(1,-1,1) , +1) -1 
1(1,-1,-1) , +1) -1 


1(2,1,0) , +1) -1 
1(2,-1,0) , +1) 1 
1(1,2,0) , +1) -1 
1(1,-2,0), +1) 1 



TABLE XIII: The momentum-space structure of Tg*", Lz = 2 sources and sinks for |np=6. These 



are shown graphicahy in figs. 14 



n ■ 



=6 



('=2) 



1(2,1,1) 
1(2,-1,1) 
1(1,2,1) 
1(1,1,2) 
1(1,-1,2) 
1(1,-2,1) 



+1) 
+1) 
+1) 
+1) 
+1) 
+1) 



1 

3 

_i 
3 

1 
3 
1 

6 
_1 
6 

_1 
3 



1(2,1,-1) 
1(2,-1,-1) 
1(1,2,-1) 
1(1,1,-2) 
1(1,-1,-2) 
1(1,-2,-1) 



+1) 
+1) 
+1) 
+1) 
+1) 
+1) 



n 



=6 



1(2,1,1) 
1(2,-1,1) 
1(1,2,1) 
1(1,1,2) 
1(1,-1,2) 
1(1,-2,1) 



+1) 
+1) 
+1) 
+1) 
+1) 
+1) 



]_ 

1 

]_ 

V2 
3 

_Vl 

3 

1 

6V2 



(<=4) 



1(2,1,-1) 
1(2,-1,-1) 
1(1,2,-1) 
1(1,1,-2) 
1(1,-1,-2) 
1(1,-2,-1) 



+1) 
+1) 
+1) 
+1) 
+1) 
+1) 



]_ 

6v^ 
1 

6V2 
]_ 

3 

_V1 

3 

1 

QV2 



that of a 2 X 2 matrix, which has solutions 
cot52 , cot54 Ao(l;g') , 7^4,0 (l;g') 10^6,0 (1;?^) 



±- 



Il7r3/2g5 Ilv/l37r3/2g7 

10^4,0 (l;g2) , 20^6,0 (l;g') 



777r3/2g5 Ilv/l37r3/2g7 



+ 4 



( 20v/3Z4,o (1; 

1 777r3/2g5 



cot52 — coiSi 



40^/ 1^^6,0 (l;g'^ 



21 1/2 



Il7r3/2g 



«7 



In the Hmit that land 4^ « tan(52 the energy-eigenvalues are the solutions to 



g^cot(52 



277 



TT 



3/2 



g%o(i;g^)-^Ao(i;g^)) = (f)'j7^'^^(^") 



(29) 



(30) 



where is shown as a function of in fig. 12 The T2 irrep first appears in the |n 



shell, as can be seen in fig. 12 The graphicaJrepresentations of the sources and sinks that 
generate the T2 irrep for the low-lying |n|2-shells are shown in fig. 13 and fig. 14, and the 
momentum-space structures for = 2 and = 1 are given explicitly in table XII and 



table |XV[ respectively. The structures of the = 3 sources and sinks are related to those 
with Lz = 1 by complex conjugation of the coefficients. 
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r+ : |Hp = 2 : L, = 2 




: |rtp = 3 : L- = 2 




FIG. 13: The momentum-space representations (left) and position-space representations (right) of 
two-body relative states in the representation with Lz = 2 for |np=2, 3, and 5 shells. 
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(1=2) : \nf = 6 : L, = 2 




(1=4) : \n\'=6 : L- = 2 




FIG. 14: The momentum-space representations (left) and position-space representations (right) of 
two-body relative states in the T2 representation with = 2 for the |np=6-shells. 



TABLE XIV: The momentum-space structure of Tg*", = 1 sources and sinks for |np=2, 3, and 
5. 



|np=2 


|n|2=3 


|np=5 


1(1,0,1), +1) -\ 
1(1,0,-1) , +1) \ 
1(0,1,1) , +1) -f 
1(0,1,-1) , +1) I 


1(1,1,1), +1) -\ 
1(1,1,-1), +1) \ 

1(1,-1,1), +1) 1 
1(1,-1,-1) , +1) 


1(2,0,1) , +1) 
1(1,0,2) , +1) 
1(0,2,1) , +1) 
1(0,1,2) , +1) 


1(2,0,-1) , +1) ^ 
1(1,0,-2) , +1) ^ 
1(0,2,-1) , +1) ^ 
1(0,1,-2) , +1) ^ 
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TABLE XV: The momentum-space structure of T2 , 



1 sources and sinks for |np=6. 



1(2,1,1) , +1) 
1(2,-1,1), +1) 
1(1,2,1) , +1) 
1(1,1,2) , +1) 
1(1,-1,2), +1) 
1(1,-2,1) , +1) 



n 



=6 



_2_ 

(T 

l_2i 

2 3_ 

2 , i 
3^ 

10 



4 I 3* 

5 5 



10 

■1«./^ 
2 



a=2) 

1(2,1,-1) , +1) 
1(2,-1,-1), +1) 
1(1,2,-1) , +1) -l^f 
1(1,1,-2), +1) 

1(1,-1,-2), +1) -i 



:i,-2,-i) , +1) 



—I 

To 



6''V 2 



1(2,1,1) , +1) 
1(2,-1,1) , +1) 
1(1,2,1) , +1) 
1(1,1,2) , +1) 
1(1,-1,2) , +1) 
1(1,-2,1) , +1) 



|np=6(;=4) 



-2i 



V65 
3 _ Si 

_ 4 3_ 

V65 

8 I 3i 
3"*" 4 

1 I Hi 
4^^ 12 

12 i 



65 



2+ 



12 
'65 



1(2,1,-1) 
1(2,-1,-1) 
1(1,2,-1) 



+1) 
.+1) 
+1) 



(1,1,-2) , +1) 
:i,-l,-2) , +1) 

:i,-2,-i) , +1) 



-2i 



V65 

3_ 81 

4 3_ 

V65 

8 I 3i 
3~'~ 4 

\/65 

1 I Hi 
_ 4^ 12 

_ 12 4 



2+ 



12 

65 



The I = 2 phase-shift was calculated from the energies of states in both the and 
irreps in recent work by Dudek et al [20]. Two states in each irrep were calculated below 
the 27r — >■ 47r inelastic threshold at the pion mass of the calculation. The contamination in 
the extraction of 82 from the higher partial-waves was estimated to be small. 



B. Negative Parity Systems 

The analysis of the odd-parity energy-levels, and their associated sources and sinks, par- 
allels that of the even-parity states. There are five negative parity irreps of the cubic 
group, , A2 , E~, Tf, and T2 with dimensions 1, 1,2,3, and 3 respectively. The energy- 
eigenvalues, sources and sinks for the negative-parity states are presented in the following 



IIIBl 




IIIB2 




IIIB3 




IIIB4 


and 


IIIB5 



As discussed previously, the A-^ irrep 



first appears relatively high in the spectrum, in the |np = 14 shell, and is sensitive to the 
/ = 9 and higher partial-waves. 



1. A2 Representation 

The energy-eigenvalues of states transforming in the A2 irrep [L^ = 2) receive contribu- 
tions only from interactions in the / = 3 partial-wave for / < 6, as presented in table [1} The 
A2 state associated with the F^^^^ block is (in the \l,m) basis) 

IA",2;3;1) = -^|3,2) - ^|3,-2) , 

and the solutions to eq. ^ from this irrep result from 



2n^ 






/ 7r3/2 






L . 


/ 7r3/2 



<l'-otSs = ( ^ ) ;i7, ( <l'2o,o (1; e) - f/Z,, (1; g-) + ^Z^, (l; g^) 



i (^~^) ' 



(31) 
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FIG. 15: The function 'Vj^^j, as defined in eq. (31), as a function of . The vertical dashed 
hnes denote the position of the poles of the function corresponding to the non-interacting energy- 
eigenvalues. 



TABLE XVI: The momentum-space structure of A2 sources and sinks. These are shown graphi- 



cally in fig. 16 




1(1,1,1),-!) \ 

:i,i,-i),-i) -\ 

:i,-i,i),-i) -\ 

1(1,-1,-1) ,-1) \ 



where the function Xj^ (q ) is shown in fig. 15 The graphical representations of the sources 



and sinks that generate the A2 irrep in the low-lying |n|2-shells (jnp = 3 and |n|2 = 6) are 



shown in fig. 16 , and the momentum-space structures are given explicitly in table XVI 

The A2 irrep first appears in the jn^ = 3-shell and / = 3 is the lowest contributing partial- 
wave. LQCD calculations of correlation functions from sources and sinks transforming as 
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FIG. 16: The momentum-space representations (left) and position-space representations (right) of 
two-body relative states in the A2 representation for the lowest-lying |np-shells. 

A2 will provide determinations of ^3 with contamination from partial-waves with / > 7, 
i.e. the energy of the A2 states receive contributions from I = 3, 7, .... This is in contrast 
to states in the T2 irrep, which will be considered subsequently, whose energy-eigenvalues 
receive contributions from partial- waves with / = 3, 5, .... This suggests that the A2 irrep is 
optimal for determining ^3. 

2. E~ Representation 

The energy-eigenvalues of E~ states receive contributions only from interactions in the 
I = 5 partial-wave for / < 6, as presented in table [l| As the E~ irrep is two-dimensional, 
the contribution to the determinant in eq. ^ results from a 2 x 2 matrix for / < 6, which 
collapses down to a one-dimensional factor as the = and = 2 states are degenerate. 
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FIG. 17: The function as defined in eq. (32), as a function of cp'. The vertical dashed 
hues denote the position of the poles of the function corresponding to the non-interacting energy- 
eigenvalues. 



state associated with the F^^^^ block is 



The E- L, 

|E-,0;5;1) 

The solution to eq. ^ from the E~ irrep results from 



27rV' 1 



L 



TT 



3/2 



-10^ A ~2N 6g'^Z4,o(l;g^) , 32q^Ze.o {1; g^) 672^^8^(1^ 1152 v^Sio,o(l;g') 

Q ^0,0 [i-',C[ ) — \ . ^ r-r _ r-^ H 



13 



17V13 



247V17 



4199 



2nV' 1 



L 



7,3/2 
;;2 



(32) 



n 



The graphical representations of the source 

2 



6-shell are shown in fig. 18, and the 



where the function X^{q ) is shown in fig. 17 
and sink that generate the E~ irrep in the 
momentum-space structure is given explicitly in table XVII 

The E~ irrep first appears in the |np = 6-shell and Z = 5 is the lowest contributing partial- 
wave. LQCD calculations of correlation functions from sources and sinks transforming as 
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TABLE XVII: The momentum-space structure of E , Lz = 0, 2 sources and sinks. The Lz = 



case is shown graphically in fig. 18 




E~ will provide determinations of ^5 with contamination from partial-waves with / > 7, i.e. 
the energy of the E~ states receive contributions from / = 5, 7, .... The LQCD calculations 
will need to be performed in relatively large volumes, as we discuss later, in order for the 
Inp = 6 shell to lie below the inelastic threshold. 



3. Representation 

The energy-eigenvalues of states transforming in the irrep receive contributions from 
interactions in the / = 1,3,5,... partial-waves, as presented in table |l| As the Tf irrep is 
three-dimensional, the contribution to eq. ^ is the determinant of a 12 x 12 matrix for / < 6 
(there are two Tf 's in the decomposition of Z = 5), which collapses down to the determinant 
of a 4 X 4 matrix as the Lz = 0, Lz = 1 and Lz = 3 states are degenerate. The Lz = 
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states associated with the F^^^\ F^r^r^\ and F^^^'' blocks are 



|Tf,0;l;l) 
|Tr,0;3;l) 
|Tf,0;5;l) 

|Tr,0;5;2) 



|i,o) 

|3,0) 
|5,0) 



|5,4) + |5,-4) 



With these four basis states, the contribution to eq. rt2|) becomes 



det 



/cot5i \ 

cot(53 

cot^s 

y cot5^J 



/p(FF,Tf) -^(FV,T-) -^(FV,T~) -^(FV,T-)\ 
^ 1;1 ^ 1;3 ^ l;5i ^ 1;52 



-p(FV,T-) p(Fy,T-) -^(FV,T-) -^(FV,T-) 

-^3;1 3;3 3;5i 3;52 

p(FV,T-) p(Fy,Tr) -^{FV,T-) -^{FV,T-) 

^ 5i;l -^5i;3 5i;5i 5i,52 

(FV.Tf) ^(Fy,T-) -^{FV,T-) -^{FV,T-) 



52;l 



52 ;3 



52,5i 



52,52 



(33) 
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where 



-^1,1 
^1,3 



^1,52 

^ 3,3 



3,5i 



F. 



3,52 

iFV,T{ 
5i,5i 

(FV,Tj- 
5i,52 

(FV,Tj- 

52,52 



7^3/2g 



2l7r3/2g5 



337r3/2g5 



7]-3/2g7 



7]-3/2g5 
I. ;^2^ 



7]-3/2g7 

I. ;^2^ 



Zo,o(l;g') , 100^6,0 ^6^4,0 



;t2\ 



7r3/2g 



33v^7r3/2g7 Il7r3/2g5 



60^4,0 (1; g^) ^ '^V 143^6,0 (1; g' ^ 



56^/4^^8,0 (1;^ 



13V777r3/2g5 



37r3/2g 



o7 



137r3/2g9 



- + - 



^8,0 



;^2> 



137r3/2g5 



7r3/2g7 



397r3/2g 



o9 



Zo,o(l;9^) 6^4,0 80^6,0 490^8,0 (1; g^) 756v/2lZio,o (1; g^) 

7r3/2g 137r3/2g5 5l^n^/^q7 247v^7r3/2g9 41997r3/2gii 

6y^Z4,o (1; e) 8^^6,0 (1; g^) 154^^8,0 (1; g^) 2772^1^10,0 (1; g^) 



137r3/2g5 



+ 



177r3/2g 



a7 



74l7r3/2g 



n9 



41997r3/2g 



nil 



Zo,o (1; g2) 6^4,0 (1; g^) _^ 32^6,0 (l;g^) ^ 14^^8,0 (1; g^) 84v^Zio,o (1; g') 



7r3/2g 



137r3/2g 



nf5 



17V137r3/2g7 



2477r3/2g 



;^9 



3237r3/2gii 
(34) 



In the hmit of vanishing interactions in the / = 5 partial- wave, eq. (33) collapses down 
to the determinant of a 2 x 2 matrix, which has solutions 



;;2^ 



cot^i cot53 2:0,0 (l;g^) 32:4,0 (l;g^) 50^6,0 (l;g' 



7r3/2g 



Il7r3/2g5 33v/l37r3/2g7 



±- 



6^4,0 (l;g2) , 100^6,0 (l;g2) 



Il7r3/2g5 33v^7r3/2g7 



+ COt^i — COt5; 



64^4,0 (l;g 
2l7r3gio 



■ (35) 



In the situation where tan^s << tan^i, eq. (35) can be perturbatively expanded to give the 
/ = 1 dominant solution 



g^cot^i 



2tt 

T 



1 



3/2 



TC 



g%o(i;g') 



27r 

T 



TT' 



3/2 



(36) 



where the function Xrp^ (g^) is shown in fig. 19 The graphical representations of the source 



and sink that generate the irrep for the lowest-lying |n| -shell are shown in fig. 20 



fig. 21 and fig. 22, and the momentum-space structure is given explicitly in table XVIII and 
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T-f : |Hp = 1 : L, = 




rf : |rtp = 2 : L- = 




Jf : |Hp = 3 : L- = 




FIG. 20: The momentum-space representations (left) and position-space representations (right) of 
two-body relative states in the Tj~ representation with Lz = for the |np = 1, 2, 3-shells. 
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T-f : |«p =4 : L, = 




Tf (1=1) : =5 : L - = 




rf (1=3) : =5 : L- = 




FIG. 21: The momentum-space representations (left) and position-space representations (right) of 
two-body relative states in the Tj~ representation with Lz = for the |np = 4, 5-shells. 
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Jf (1=1) : \nf = 6 : L - = 




Tf (1=3) : \n\' =6 : L - = 




FIG. 22: The momentum-space representations (left) and position-space representations (right) of 
two-body relative states in the Tj~ representation with = for the |np = 6-shells. 



TABLE XVIII: The momentum-space structure of , = sources and sinks for |np=l-4. 



These are shown graphically in figs. 20 and 21 



|n|2=l 


|n|2=2 


|n|2=3 


|n|2=4 


1(0,0,1) ,-1) 1 


1(1,0,1), -1) ^ 
1(1,0,-1) ,-1) -1 

1(0,1,1), -1) I 
1(0,1,-1) ,-1) -1 


1(1,1,1), -1) -h 

1(1,1,-1) ,-1) 5 
1(1,-1,1) ,-1) -1 

1(1,-1,-1) ,-1) 2 


1(0,0,2) ,-1) 1 
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TABLE XIX: The momentum-space structure of T-^ , 



sources and sinks for |nP=5 and 6. 



These are shown graphically in figs. 21 and 22 



|n|2 




1) 


|n|2: 




3) 


|n|2 = 


-6/; 


1) 


|n|2 = 


-6/; 


3) 
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1 


1 
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1 V ' ' / ' 
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1 

2^3 














1(2,1,-1) 


, -1) 


1 

2v^ 


1(2,1,-1) 

IV'' / 


, -1) 


1 

2^3 


1(2,0,1) 


-1) 


1 

2V5 


1(2,0,1) 


-1) 

1 


1 


1(2,-1,1) 

IV' ' / 


, -1) 


1 

2v^ 


1(2,-1,1) 

IV' 5 / 


, -1) 


1 


1(2,0,-1) 


,-1) 


1 


1(2,0,-1) 


,-1) 


1 


1(2,-1,-1) 


,-1) 


1 

2v^ 


1(2,-1,-1) 


,-1) 


1 

2V3 


1(1,0,2) 


-1) 


1 

V5 


1(1,0,2) 


-1) 


1 

2v^ 


1(1,2,1) , 


-1) 


1 

2^6 


1(1,2,1) , 


-1) 


1 

2V3 


1(1,0,-2) 


, -1) 


1_ 

V5 


1(1,0,-2) 


, -1) 


1_ 

2V5 


1(1,2,-1) 


, -1) 


1 

2v^ 


1(1,2,-1) 


, -1) 


]_ 

2^3 


1(0,2,1) 


-1) 


1 


1(0,2,1) 


-1) 


1 


1(1,1,2) , 


-1) 


1 

V6 


1(1,1,2) , 


-1) 


1 

2^3 


1(0,2,-1) 


,-l> 


1 


1(0,2,-1) 


,-1) 


1 

V5 


1(1,1,-2) 


,-1) 


1 


1(1,1,-2) 


,-1) 


1 

2^3 
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-1) 


1 
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-1) 
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V6 
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1 

2^3 
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,-1) 


1 

V5 
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,-1) 


1 

2v^ 


1(1,-1,-2) 
1(1,-2,1) 
1(1,-2,-1) 


,-1) 
,-1) 
,-1) 


1 

1 

2^6 
1 

2v^ 
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1(1,-2,1) 
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,-1) 
,-1) 
,-1) 


1 

2^/3 
1 

2V3 
1 

2V3 



TABLE XX: The momentum-space structure of , = 1 sources and sinks for |np=l-4. 



|n|2=l 


|n|2=2 


|n|2=3 


|n|2=4 


1(1,0,0) ,-1) ^ 
1(0,1,0) ,-1) ^ 


1(1,1,0), -1) ^ 
1(1,0,1), -1) 1 + 1 
1(1,0,-1) ,-1) 1 + 1 

1(1,-1,0) ,-i> I 
1(0,1,1), -1) -l + i 
1(0,1,-1) ,-1) -l + i 


1(1,1,1), -1) 2 
1(1,1,-1) ,-1) I 
1(1,-1,1) ,-1) --2 
1(1,-1,-1) ,-1) 


1(2,0,0) ,-1) ^ 
1(0,2,0) ,-1) ^ 



table |XX[ The structures of the = 3 sources and sinks are related to those with = 1 
by complex conjugation of the coefficients. 

The T]~ irrep first appears in the |np = 1-shell and / = 1 is the lowest contributing partial- 
wave. LQCD calculations of correlation functions from sources and sinks transforming as 
will provide determinations of 61 with contamination from partial-waves with / > 3. LQCD 
calculations of the phase-shift in this partial-wave are presently being performed, and the 
p-resonance is beginning to be mapped out, e.g. Ref. ^8] . 



4- T2 Representation 

The energy-eigenvalues of states transforming in the Tg" irrep receive contributions from 
interactions in the / = 3, 5, ... partial-waves, as presented in tablejlj As the T2 irrep is three- 
dimensional, the contribution to the determinant in eq. ^ results from a 6 x 6 matrix for 
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TABLE XXI: The momentum-space structure of T-^ , 



1 sources and sinks for |np=5 and 6. 



|np=5(;. 



1(2,1,0) 
1(2,0,1) 
1(2,0,-1) 
1(2,-1,0) 
1(1,2,0) 
1(1,0,2) 
1(1,0,-2) 
1(1,-2,0) 
1(0,2,1) 
1(0,2,-1) 
1(0,1,2) 
1(0,1,-2) 



1 

2V2 
2 _i 

5 5_ 

V2 

2 _1 

5 5_ 

V2 

3__2i 
10 5 

V2 
2 ,31 
5"*" 10 

1^- 

1 7^ 

5 10 

V2 

1 i_ 

5 10 

_ i 

1 ,2i 
5"*" 5 

V2 
1 I 2i 
5"*" 5 

V2 

j_ , 1 
10^5 

j_ , 1 
10^5 



V2 



n 



^(^=3) 



(2,1,0) , 
(2,0,1) , 
(2,0,-1) 
(2,-1,0) 
(1,2,0) , 
(1,0,2) , 



(1,-2,0) 
(0,2,1) , 
(0,2,-1) 
(0,1,2) , 
(0,1,-2) 



-1) 
-1) 
-1) 
-1) 
-1) 
-1) 



(1,0,-2) , -1) 



-1) 
-1) 
-1) 
-1) 
-1) 



1 

J_-i_i 
10^5 

V2 

, i 
10^5 



2i 



10 



V2 
2 I 3i 

5"*" 10 

V2 
1 I 2i 
5"*" 5 

1 I 2i 
5"*" 5 

_ i 
2V2 



^V2. 

_5 10 

V2 

2_i 

5 5_ 

5 5. 

V2 



n 



=6 



('=!) 



1(2,1,1) ,- 

1(2,1,-1) ,- 
1(2,-1,1) ,- 

1(2,-1,-1) , 
1(1,2,1) ,- 
1(1,2,-1) ,- 
1(1,1,2) ,- 
1(1,1,-2) ,- 
1(1,-1,2) ,- 
1(1,-1,-2) , 
1(1,-2,1) ,- 
1(1,-2,-1) , 



1) 

■1) 

■1) 

-1) 

1) 

•1) 

1) 

■1) 

■1) 

-1) 

■1) 

-1) 



/ 5 
V 3 
4 

4 _ 

15 



V15 
15 



15 



_±_ _L _i_ 

15 10 

A 4_ _L 
15 10 



_2_ 
15 



10 



2^ 
15 



10 



n 



^(^=3) 



1(2,1,1) , 


-1) 


1(2,1,-1) 


-1) 


1(2,-1,1) 


-1) 


1(2,-1,-1) 


,-1) 


1(1,2,1) , 


-1) 


1(1,2,-1) 


-1) 


1(1,1,2) , 


-1) 


1(1,1,-2) 


-1) 


1(1,-1,2) 


-1) 


1(1,-1,-2) 


,-1) 


1(1,-2,1) 


-1) 


1(1,-2,-1) 


,-1) 



2^3 
i 

2^3 
]_ 

2^3 
]_ 

2^3 
_ i 

2^3 
_ i 

2^3 
1 

2v^ 
1 

2v/3 
_ i 

2^3 

_ i 

2v^ 
1 

2y/Z 
1 

2^3 



Z < 6, which collapses down to the determinant of a 2 x 2 matrix as the = 1, = 2 and 

L2 = 3 states are degenerate. The = 2 states associated with the F^^^^ and -^5^5^^^ 

blocks are 



^13,2) + i=|3,-2) 



^15,2) + i=|5,-2) 



K,2;3;l) 
|T2-,2;6;1> 

in terms of which, the Tg" contribution to eq. ([2|) becomes 
det 



cot (^3 
cot(55 



-piFV,T^) -^{FV,T^)^ 

3"3 35 

-^{FV,T-) -^(fV,T-) 

^ 5;3 ^ 5;5 



where 

^{FV,T-) 



F 

F. 
F 



3;3 

iFV,T^) 
3;5 

iFV,Ti) 
5;5 



2o,o{l;e ) 2^4,0 (l;g2) 60^6,0 (l;g 
20^4,0 (l;g 



;^2\ 



;;2\ 



Il7r3/2g5 11^7r3/2g7 

14^6,0 (l;g2) ^ 112^8,0 (l;g') 



(37) 



13v^7r3/2g5 v^Vr3/2g7 13v^7r3/2g9 

Zo,o (1; g') ^ 4^4,0 (1; q') 80Ze,o (1; g') 280^8,0 (1; g^) 432y2lZio,o (1; g') 



137r3/2g5 17v^7r3/2g7 247v^7r3/2g9 



41997r3/2fi 
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TABLE XXII: The momentum-space structure of , = 2 sources and sinks. These are shown 



graphicahy in figs. 24 and 25 



n 



1(1,0,1) ,-!> 

1(1,0,-1) ,-i: 

1(0,1,1) ,-1) 

1(0,1,-1) ,-i; 



|np=5(;=3) 



1(2,0,1) 
1(2,0,-1) 
1(1,0,2) 
1(1,0,-2) 
1(0,2,1) 
1(0,2,-1) 
1(0,1,2) 
1(0,1,-2) 



-1) 
-1) 
-1) 
-1) 
-1) 
-1) 
-1) 
-1) 



]_ 

1 

V5 
]_ 

1 

2v^ 
1 

1_ 

V5 
1 

]_ 

2v^ 



1(2,0,1) 
1(2,0,-1) 
1(1,0,2) 
1(1,0,-2) 
1(0,2,1) 
1(0,2,-1) 
1(0,1,2) 
1(0,1,-2) 



'('=5) 

2^5 
1 

'2^5 
1 

1 

V5 
1 

'2^5 
1 

2v^ 
1 

V5 
1 



-1 

-i; 
-1 

-i; 
-1 

-i; 
-1 



1(2,1,1) ,-1) 

1(2,1,-1) ,-i; 
1(2,-1,1) ,-i; 

1(2,-1,-1) ,-1 
1(1,2,1) ,-1) 

1(1,2,-1) ,-i; 
1(1,-2,1) ,-i; 

1(1,-2,-1) ,-1 



]_ 

2^2 

1 

]_ 

1 

2V2 
1 

2V2 

]_ 

2V2 
1 

2V2 

]_ 

2V2 



The solutions to this equation result from 
cot53 , cot55 ^0,0 (1; q^) 9^4,0 (1; <?') , 9502^,0 (1; g') , 140^8,0 (1; ^') , 216v^Zio,o (1; g') 



3/2 



g 1437r3/2g5 187v^7r3/2g7 247v^7r3/2g9 41997r3/2g 



nil 



±- 



432x/21Zio,o (1; g') , 70^4,0 (l;g') , 140^6,0 (1; g') 280^8,0 (1; g'^ 



41997r3/2gii 



+ 



+ 



1437r3/2g5 187v^7r3/2g7 247v^7r3/2g9 



+ COtfJs — cot^s 



+ 



4 /20Z4,o(l;g~2) , 14Z6,o(l;g2) 112^8,0 (1;^ 



ISVllg 



nf5 



+ 



^l43g' 



n7 



1/2 



(38) 



In the limit of vanishing interactions in the / = 5 partial- wave, eq. (37) collapses down to 

7 



g''cot53 



2n 

T 

271 

T 



1 



IT 



3/2 



3/2 



TC' 



g%o(i;g')- 



^T. if) 



2^2^4,0 (l;g2) 60^6,0 (l;g'^ 



11 



iiyis 



(39) 



where the function Xj,^ (g^) is shown in fig. 23 The graphical representations of the source 
and sink that generate the T2 irrep for the lowest-lying |np-shell are shown in fig. 



24 



fig. 25 , and the momentum-space structure is given explicitly in tables XXII and table 



and 



XXIII 



The structures of the = 3 sources and sinks are related to those with = 1 by complex 
conjugation of the coefficients. 

The T2 irrep first appears in the |np = 2-shell and / = 3 is the lowest contributing partial- 
wave. LQCD calculations of correlation functions from sources and sinks transforming as 
T2 will provide determinations of 63 with contamination from partial-waves with / > 5. 



5. A-^ Representation 

The lowest-lying state transforming in the irrep is in the |np = 14-shell. The energy- 
eigenvalues are sensitive to interactions in odd partial- waves with / > 9, and the energy- 
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FIG. 23: The function Af^^' ^ defined in eq. (39), as a function of q . The vertical dashed 



hnes denote the position of the poles of the function corresponding to the non-interacting energy- 
eigenvalues. 



splitting in the large volume limit is dominated by the / = 9 partial-wave. Using the methods 
of the previous section to isolate the state and determine the appropriate energy-eigenvalue 
equation is tedious as i^g.g^^'' is a 19 x 19 matrix. Using the following spherical basis state, 

l^r,0;9; 1) = ^y|(|9,8) - |9,-8)) - " |9, -4)) , (40) 
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77 : |Hp = 2 : L, = 2 




n (l=3>) : \n\'=5 : L- = 2 




T2 (1=5) : =5 : L- = 2 




FIG. 24: The momentum-space representations (left) and position-space representations (right) of 
two-body relative states in the T2 representation with = 2 for select |np-shells. 
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77 : |«p = 6 : L, = 2 




FIG. 25: The momentum space representations (left) and position-space representations (right) of 
two-body relative states in the T2 representation with Lz = 2 for the |np = 6-shell. 



TABLE XXIII: The momentum-space structure of T2 , = 1 sources and sinks. 
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2=2 






|n|2 




=3) 


|n|2 
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V2 
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3 1 2i 
lO''" 5 

V2 
2 3i 
5 10 


1(2,-1,1) , 


-1) 


1 
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i 
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\ + 


i 
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5 5 

V2 
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-1) 
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,V2. 
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-1) 


i 
4 
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-1) 


i 
2 
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1 2i 
5 5 
V2 
i 

1, j_ 


1(1,0,-2) 
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V2 

i 
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_i 1 i 
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-1) 


111 
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-5e+1 1 



-1e+12 




FIG. 26: The function X^^, as defined in eq. (41), as a function of g^. The vertical dashed hne 
denotes the position of the first pole of the function corresponding to the non-interacting energy- 
eigenvalue. 



the eigenvalue-equation for the interaction in the / = 9 partial-wave is 

T 



q cot ^9 



TT 



3/2 



32v/l3gi2z6,o(l;g2 
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115 



56Vl7q'^Zs,o{l-,_fl ^ 1568v/7g^Zio,o(l;g^) ^ 308g6Zi2,o(l; g') 
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3565v^ 
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T 



19 



TT 



3/2 



(41) 
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n 



where the function (g^) is shown in fig. 
and sink that generate the irrep for the 
momentum-space structure is given explicitly in table |XXIV[ It is interesting to note that 
these odd-parity singlet states require the integers comprising the integer triplet to differ 



The graphical representations of the source 
P = 14-shell are shown in fig. 27, and the 
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TABLE XXIV: The momentum-space structure of , Lz = sources and sinks. These 
are shown graphically in fig. [27j The coefficients of the state vectors are of the form Cj ~ 




from each other. The first |np-shell for which this is possible has |np = 14, and the next 
has |np = 21. Given the first appearance of this irrep is high in the spectrum, a LQCD 
calculation of the / = 9 phase-shift will require enormous lattice volumes in order for the 
state to lie below inelastic thresholds. Thus, this calculation cannot be expected to be 
performed in the near future. 
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TABLE XXV: The |n|2-shell of the lowest-lying energy-eigenstate transforming as r^*\ and the 
angular- momentum of the dominant interaction in the large- volume limit. 



|n|2 
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At 
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rr 
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E+ 
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T+ 
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T2 


3 


3 
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5 
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IV. DISCUSSION 

A. Strategy for Extracting Phase-Shifts from Lattice QCD 

In Lattice QCD calculations, sources and sinks with the quantum numbers of the hadronic 
states of interest generate correlation functions, which in general are sums of exponentials 
with arguments that depend upon the energy of the eigenstates in the lattice volume. One 
path for LQCD calculations to follow is to form sources and/or sinks that transform as 
irreps of the cubic group from the eigenstates of linear-momentum (generated by the Fourier- 
transform of single-hadron objects). Clearly, such sources are not the energy-eigenstates in 
the lattice-volume due to the interactions between the particles, and as such these sources 
and sinks will couple, in principle, to all states in the lattice-volume with the appropriate 
quantum numbers. At large times, the correlation function will depend exponentially upon 
the energy of the lowest eigenstate. For instance, the correlation function resulting from 
an A\ source constructed from eigenstates with |np ^ will depend upon the energy of 
the lowest-lying At irrep at large times. Of course, the overlap onto the lowest-lying At 
irrep may be small, in which case the ground-state dominates only after a large number of 
time-slices. It is then clear that the optimal way to extract the bi for Z < 6 is to determine 
the lowest energy-eigenvalue of the cubic irrep that is dominated by each 81. These irreps 



are shown in table XXV, along with the |n| -shell of the lowest-lying energy-eigenstate that 



contributes to the corresponding correlation functions. Table XXV shows that with just the 
lowest two |np-shells, |np = 0, 1, the phase-shifts in the lowest three partial-waves, 5o,i,2; 
can be determined. In order to determine the phase-shifts in all partial- waves with / < 6, 
correlation functions must be formed for states that have ground-states in the shells up to 
|np = 6. Determining the phase-shift for / > 7 can be seen to be decidedly more difficult 
than for / < 6 as there is only one further irrep of the cubic group, the A~{ which first occurs 
in the |np = 14 and is dominated by the interactions in the / = 9 partial-wave in the large 
volume limit. 

For shells in which there are multiple occurrences of a given T*^*-*, the partial- 
diagonalization of the states in the infinite-volume limit in terms of the angular-momentum 
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of the interactions is possible. However, sources cannot be constructed to isolate these states 
due to interactions, and in general, closely spaced-states will be encountered in the spectrum. 



B. Expectations for the vrvr Energy-Eigenvalues 

In order to estimate the computational resources required to extract the vrvr phase-shifts in 
higher partial-waves, the experimentally determined (and parameterized) phase-shifts can be 
used to determine the energy-eigenvalues for a range of lattice volumes. The vrvr phase-shifts 
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FIG. 28: The tttt phase-shifts, Sf , as a function of y/s < Amj^ for / = 0, 1, 2, 3, as parameterized in 
Ref. 1571. 



for / = 0, 1, 2, and 3 partial-waves extracted from experimental data and parameterized with 
functions that satisfy unitarity and analyticity, and specifically incorporate any lowest-lying 
resonances in the channel ^ are shown in fig. 
eters describing each partial-wave provided in Ref. 



28 



The central values of the param- 
are used in the analysis, but, as 
only estimates of the energy-eigenvalues are being explored, a systematic propagation of the 
uncertainties has not been performed. Further, we assume isospin symmetry in our analy- 
sis. As Liischer's formalism is valid below inelastic thresholds, only the phase-shifts in the 
kinematic regime ^/s < AitLt^ (at the physical pion mass) are considered. 



Figure 29 shows the energy-eigenvalues associated with two non-interacting pions with 
vanishing total momentum in the lattice volume (also shown in fig. 1 of Ref. [T]). At the 
physical pion mass, LQCD calculations in volumes with 6 fm are highly desirable in 
order to suppress the exponential corrections that are not included in the formalism of 
Liischer [55]. The energy-shift between the non-interacting state and the interacting 1 = 
Ai and 1 = 2 states in the |np = shell are shown in fig. 30 
the eigenstates in the |np = 1,2,3,4 shells are shown in figs. 31 34, respectively, 
energy-shifts for the |np 



The energy-shifts for 
The 



A]^ states and the |n| = 1 state can also be found in fig. 7 
of Ref. [1] . The energy-shifts of the states due to the s-wave and p-wave interactions are of 
comparable size. As the s-wave interactions are currently being calculated in volumes with 
L ~ 3.5 fm, we do not anticipate significant difficulty in performing these calculations at the 



^ The real-part of the inverse scattering amplitude, cotSj , is expanded as a power-series in the function 



w{s) 



= , where Si is the energy above which inelastic processes cannot be neglected. 
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L (fm) 

FIG. 29: The energy-eigenvalues associated with two non- interacting pions as a function of the 
spatial-extent of the lattice volume, L. The spectrum results from the momentum of each pion 
being restricted to q = for all possible triplets of integers, n, due to the periodic boundary 
conditions imposed on the quark and gluon fields. The maximum value shown on the vertical-axis 
corresponds to the inelastic threshold ^/s = Am-,^. For reference, a spatial-extent of L = 6 fm 
corresponds to m-,jL ~ 4.2. 
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L (fm) L (fm) 

FIG. 30: The expected vrvr energy-shifts in the |np = shell due to strong interactions. The 
left panel shows the shift in the / = irrep (dominated by 5q), while the right panel shows 
the shift in the 1 = 2 Af irrep (dominated by 6q). Both the L-axis and the AiiJ-axis are scaled 
logarithmically (log;^g). 

physical pion mass in lattices with L> 6 fm. In contrast, the energy-shifts of states due to the 
d-wave (/ = 2) and f-wave (/ = 3) interactions are more than an order of magnitude smaller 
than those of the Af irrep. Significantly more computational resources will be required to 
extract the phase-shifts beyond the s-wave and p-wave. It is difficult to make estimates for 
the energy-shifts due to interactions beyond the f-wave as the experimental measurements 
of these phase-shifts have large uncertainties or are absent. Given the results obtained for 
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FIG. 31: The expected tttt energy-shifts in the |np = 1 shell due to strong interactions. The 
left panel shows the shift in the / = Af , irreps (dominated by 6q and 82, respectively), the 
center panel shows the shift in the 1 = 1 irrep (dominated by 61) and the right panel shows 
the shift in the / = 2 A+, E+ irreps (dominated by 6q and (5|, respectively). Both the L-axis and 
the -axis are scaled logarithmically (logj^g). 
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FIG. 32: The expected tttt energy-shifts in the |np = 2 shell due to strong interactions. The left 
panel shows the shift in the / = Af , E'^, T2 irreps (dominated by 5q, 82 and ^2, respectively), 
the center panel shows the shift in the 1 = 1 , T2 irrep (dominated by 5\ and 6\) and the right 
panel shows the shift in the 1 = 2 A^, E^ ^ T2 irreps (dominated by 5q, 62 and respectively). 
Both the L-axis and the A£^-axis are scaled logarithmically (log^^o)- 




FIG. 33: The expected vrvr energy-shifts in the |np = 3 shell due to strong interactions. The left 
panel shows the shift in the / = Af , irreps (dominated by Sq and 82, respectively), the center 
panel shows the shift in the 1=1 , A2 irreps (dominated by 6l and 6^, respectively) and the 
right panel shows the shift in the 1 = 2 A'^ , T2 irreps (dominated by 5q and respectively). 
Both the L-axis and the Ai?-axis are scaled logarithmically (log^^g)- 



/ < 3, it is not difficult to speculate as to the size of the energy-shifts of partial-waves beyond 
/ = 3, and the associated difficulty in their extraction from LQCD calculations. 

In order to estimate the amount of mixing of higher partial waves to a given phase-shift 
from the energy-eigenvalues, it is important to understand the expected contributions from 
(all of) the partial-waves. The energy-splitting of the irrep in the |np = 1-shell from 
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FIG. 34: The expected tttt energy-shifts in the |np = 4 sheh due to strong interactions. The 
left panel shows the shift in the / = , irreps (dominated by 6q and 82, respectively), the 
center panel shows the shift in the 1 = 1 irrep (dominated by 61) and the right panel shows 
the shift in the / = 2 A+, E+ irreps (dominated by 6q and (5|, respectively). Both the L-axis and 
the -axis are scaled logarithmically (logj^g). 
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FIG. 35: The contributions to the energy-splitting of the irrep in the |np 
and S^. Both the L-axis and the A£'-axis are scaled logarithmically (logj^g). 



1-shell due to (5| 



the / = 1 phase-shift, 6l, and the / = 3 phase-shift, ^3, are shown in fig. 35, As expected, 
the contribution from ^3 is approximately two-orders of magnitude smaller than that from 
6l over the range of lattice volumes for which the analysis is applicable. Therefore, to high 
precision it is sufficient to use the perturbative expansion of the energy-splitting in terms of 
tan 61- 



C. Signal-to-Noise Issues 

There is also a signal-to-noise "problem" in the extraction of the Si for / > 1 as the signal-to- 
noise ratio degrades exponentially at large times. To demonstrate this behavior we return 
to the argument given by Lepage jSS]. Consider the correlation function resulting from a 
source that creates a 7r+7r+-state that transforms in the irrep of the cubic group, 

xit) = {9EAt)) = (0| SEAt)SU{0) |0) Ze+ e-''o^+^'* , (42) 

where SE+{t) annihilates a 7r+7r+ in the irrep at the time t. At large times this corre- 
lation function depends exponentially upon the ground-state energy which, in the absence 

of interactions, is E^^^+ ■* = 2\/ (^)^ + ml.. The variance of this correlation function is 



54 



given by 



(0| 5^+(t)44t)4^(0)5^+(0) |0) - 
J2 Cr (0| 5r(t)5r(0) |0) - 



— )• e 



-2E 



- 72 



-2E 



0.E+ ^ 



(43) 



where Cr are the Clebsch-Gordan coefficients in the expansion ® = © © 
|53l [591 EQ]. The energy, Ans, that dictates the long-time behavior of the variance 
correlation function is that of the lowest-lying irrep composed of four pions, which is the 
lowest-lying A^-irrep that has an energy of 2E^^^ ^ = 4m^ in the absence of interactions. 
Therefore, at large times, the noise-to-signal ratio behaves as 



a(t) 



^ A + 

1 gANSt 



^E+ 




(44) 



which grows exponentially at large times. 

This argument generalizes to all of the cubic irreps, and the extraction of the 61 for each 
/ > 1 suffers from a signal-to-noise problem, with an energy-scale that is approximately 



,(H= 




n r -|- 



(45) 



Obviously, in large volumes, the degradation of the signal obtained for low-lying states will 
not dramatically impact the determination of the energy-eigenvalues as the energy-scale 



behaves as A 



(|n|2) 
NS 



|np/(m7rL^). However, for present-day calculations in modest lattice 



volumes, the degradation of the signal may impact the extraction of the phase-shifts in 
higher partial- waves, and numerical exploration is required to determine its impact. 



V. CONCLUSION 

We have explored the phenomenology of Liischer's method in the extraction of the phase- 
shifts in higher partial-waves describing meson-meson scattering below inelastic thresholds 
using Lattice QCD, the formalism for which is contained in the works of Liischer. The 
lowest-lying s-wave and p-wave interactions were explored in those works, and at the time, 
Lattice QCD calculations of scattering beyond the s-wave and p-wave were in the distant 
future. However, the rapidly increasing computational resources that are being directed 
towards Lattice QCD calculations will allow for the calculation of the phase-shifts in higher 
partial-waves in both mesonic and baryonic processes in the near future. We have considered 
the low-lying spectrum of two-meson states in a finite cubic lattice-volume, and determined 
contributions from the partial-waves with I < 6 and / = 9. There are a sufficient number 
of irreps of the cubic group that will allow for the calculation of the phase-shifts, 61 for 
I < 6, and possibly I = 9. There are no irreps of the cubic group with ground-state energy- 
splittings that are dominated by interactions in the / = 7 and / = 8 partial-waves. As 
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such, there appears to be no clean way to calculate these phase-shifts from Lattice QCD 
calculations performed in cubic volumes. High precision calculations of the energy of states 
in other irreps may allow for their extraction by forming differences of energies, but this 
will require significantly more computational resources than the extraction of the phase- 
shifts in lower partial-waves. We have provided the structure of sources that will produce 
the irreps of the cubic group, in both momentum and position space, that will generate 
the relevant states in LQCD calculations. Further, we have given the explicit formula, 
and their perturbative solutions, that are required to analyze the results of such LQCD 
calculations. We recapitulate the leading contributions to the energy-eigenvalue equations 
and their solutions in the large- volume limit in table XXVI and table XXVII respectively. 

Experimental measurements of the vrvr phase-shifts are difficult, with precision currently 
at the few-percent level in the s-wave and p-wave, but much larger uncertainties are as- 
sociated with the phase-shifts in the higher partial-waves. It would appear that Lattice 
QCD calculations will be able to provide low-energy meson-meson phase-shifts in the low 
partial-waves with significantly more precision than the corresponding experimental mea- 
surements. While the contributions to energy-splittings rapidly become smaller with higher 
partial-waves, we conclude that it is presently possible to extract phase-shifts for partial- 
waves with / < 3. The implications of the recent preliminary calculations [20] of the / = 2 
phase-shift, at unphysical pion masses are very encouraging for future calculations. 



Appendix A: Block-Diagonalization of F^^^^^ 

As the number of partial-waves with non-zero phase-shifts increases, so does the complexity 
of the calculation of the energy-eigenvalues in a finite cubic volume. To illustrate the method 
for determining the energy-eigenvalues, we provide the details of the calculation when 5; 7^ 
for / < 4. As is true in all cases involving parity-conserving interactions, the analysis in 
the even-parity sector (/ = 0,2,4) decouples from that in the odd-parity sector (/ = 1,3). 
The calculations that are required for 5; 7^ for / > 4 become more complicated due to the 
dimensionality of the matrices involved, the contributions from the Zi^m for higher values 
of /, and to multiple occurrences of the same cubic irreps. This last feature means that the 
diagonalization of the finite- volume functions, Fif^\ are not dictated entirely by geometry 
due to mixing between the multiple occurrences of a given T^^\ However, such calculations 
are a straightforward extension of what follows. 



1. Odd-Parity Sector with ^1,3 / 



In the odd-parity sector with only 5i 3 7^ 0, the finite- volume corrections are encapsulated 
in which is a 10 x 10 matrix. It has block form 





^i;3 




^3;3 



(Al) 
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TABLE XXVI: A summary of the energy-eigenvalue equations for the lowest-lying state in each 
r*^*-* arising from the interaction in the dominant partial-wave. 



I [T^ 


Leading Eigenvalue Equation 


Section 


0[At] 


gcot(5o = ^0,0 (1 






IIIAl 




1 [Tf ] 




%o(l;g^) 




IIIB3 




2[E^] 
2 [T+] 


q cotd2 - (xJ ^ I 


g4Zo,o(i;g2) + fAo(i;9^)) 

;;4-7 /'i.;;2\ 4-7 /'i.;j2\\ 




III A 3 

TTT A ^ 
ill A 




3 [T^] 
3 [A,] 


=(^)'^( 


^gZo,o(l,gj 11 j 

>^o,o (1; ?~^) (1; 9^) + i/;i3^6,o (1; e)) 




IIIB4 

TTT T ") 1 

IIIB 1 




4 [r+] 


9'cot54 =(^)'^ 


in^'^ h■;^2^ 448^8,o(l;g2) 4g-226,o(l;g2) 54g424.o(l;g-2) \ 
(^(/Zo,oil,gj ^43^ + 143 j 




III A 4 




5[E-] 


,"cot<55 =(^)"^x 

/^~10^ ^^ ~2\ 6g«24,o(l;g') , 32g426,o(l;g2) 672g2^8,o(l;'?2) , Il52v^2io.o(l;g2) ^ 
^0,0(,-l-,9j 13 1 247v^ ' 4199 ) 




IIIB 2 






/^~12 7 (^ ~2\ , 6g«Z4,o(l;g2) l60v^g'3^6,o(l;g^) 40q*^8,o(l;g2) 
lg ^0,Oii,Q'j+ 17 323 19^ 

2592V2Tg2 2-^o_o(i;q;2) 19802i2,o 264v^T00T2i2,4 (l;^^) \ 
7429 ' 7429 ' 7429 ) 




TTT A 9 




9 [A^] 


cot ^9 = (f)^^ 

(^g ^0,0(1, ^ ) 23 115 345 

1568v^Q«2io,o(l;g^) 1 308g62i2,o{l;g2) , 616v^lOOTg62i2,4(l;g2) , 53248^4 Zi4_o{l;52) 

^ 3335^3 ' 2139 ' 20677 ' 10695v^ 

1664v^g2_2^g jj(l.g2) 832V46189g22ig 4(l;g2) 2206464Zi8,o{l;'?2) 28288v^35532i8,4(l;g^) \ 

3565v^ 103385v4 103385\/37 20677^25^ / 




IIIB 5 
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TABLE XXVII: The perturbative expansions of the energy-eigenvalues of the lowest-lying state in 
eachrW. 



|np [r»] 




0[^+] 


\l 2.8373 ('^^0) +6.3752 ('^^0)^ 1 + ... 


1 [rrl 

1 [E+] 


[ 1 - ^tan(5i (1 - 0.3653 tan 5i) + ...] 
^ [ 1- ^tan 52 (1 - 1.5672 tan ^2) + ...] 


2 [T+] 


47r2 
IX 


2 2X^t^^<^2(l 
[2 sX^tan^a + ...] 


4830tan(52) + ... 


3 [^2"] 


47r2 

~w 


[3 9Z^tan<53 + ...] 




5 [r+] 


IX 
47r2 

TX 


"^-lil^tan^, + ...] 
.5-3Sfetan5a + .. 




6[E-] 


4tt'^ 

IX 


[6 l2'X^tan55 + ...] 




14 [A^] 


TX 


1 ^ 4208972625 x 1 1 

[ 1^ 46118408^14^2 tan 59 + ...J 
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where the component matrices in the |/,m)-basis are 



^1:1 



-^1;3 



-,{FV) 
3;3 



3/2 



; 3o,o(l;r) diag (1,1,1) 



7]-3/2 








i 








2 




















3/2 



3 Ao(l;g') diag (1,1, 1,1, 1,1,1) 



+ 



y^3/2 g5 



— 2^4,0(1; g 



;;2> 



+ 



1 1 



7^3/2 ql 



^6,0 





/ 


3 















\ 















-7 








5 









1 










1 





A 


/Is 






11 










6 

















/Is 








1 



















5 








-7 











\ 





A 


/Is 








3 ) 












/ 


-1 








7 


/Is 





\ 













6 











-42 





5 
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-15 











7\/l5 












^15 










20 

- 




-15 










33/ 


13 




7 













-42 











6 











V 





7 













-1 / 



(A2) 



where the relevant relations between the Zi^m, that can be found in eq. (B7), have been 
used, and where diag(a, 6, ...) denotes a diagonal matrix. It is convenient to first diagonalize 



the F;^f blocks (Fi^fi 



{FV) 



is already diagonal in this basis). The block-diagonal matrix, 



is defined to have the form 



Six 








5'33 



(A3) 



and when acting on i^^^^ produces a matrix, F^^^ = S^.f[^^\sU which can be re- 
arranged into block-diagonal form where each block is associated with a r^*\ The matrices 
cos 5 and sin 5 in eq. ([3]) are invariant under this transformation. 



COs5 = diag(ci,Ci,Ci,C3,C3,C3,C3,C3,C3,C3) 



(A4) 



Explicitly, the basis is {|1, 1), |1,0), |1, -1) , |3, 3), |3,2), |3, 1), |3,0), |3, -1), |3, -2), |3, -3)}. 
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of S- in eq. (A3) are 



where Ci and C3 denote cos Si and cos ^3, respectively, and similarly for sinS. The components 

/ n n ^\ 

^0 
I 



^11 




33 



and the matrix F^^^^ is of the form 



-,(FV) 



( 






- 





1 




u 





1 








^^\ 











1;1 


^i;3 


^3;1 


3;3 







^0 






1 / 



I 



(A5) 



(A6) 



where 



^ 3:3 



7]-3/2 q 



3Zo,o(l;r) diag (1,1,1) 



1 1 



7^3/2 g5 



;;2^ 



0-1 
1 
-1 



^Zo,o(l;g') diag (1,1, 1,1, 1,1,1) 
+ ^TTT 4 Ao(l;g') ^ diag (-1,-1, -1,-6, 3, 3, 3) 



7^3/2 g5 



H ^ t:^ ^6,0(1; ? 



From the structure of F^^"* it is clear that the ordering of the F^*^ along the diagonal of 
F^^^^ is T2 , A2 and Tf, respectively, and the equations that dictate the energy-eigenvalues 
of each of the F'^*\ given in eqs. (31), (35), and (39), follow directly from eq. (A7). The 
matrix 5*- that diagonalizes F^'^^ is independent of the Zi^m functions because, with each 
relevant F^*^ occurring only once, the decomposition depends upon geometry only. 



;;2^ 



20 



33V13 



diag (-9, -9, -9, 12, 5, 5, 5) 



(A7) 



60 



2. Even- Parity Sector with (5o,2,4 7^ 



In the even-parity sector with only (5o,2,4 7^ 0, the finite- volume corrections are encapsulated 



in F_ 



(FV) 



which is a 15 X 15 matrix. It has block form 



where 



p{FV) 





/ 
( 

j^iFV) 


T^{FV) 
^Q;2 


\ 

p{FV) 




piFV) ^ 


piFV) 
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^2-2 


piFV) 
^2;4 


(A8) 




-^4;0 

\ 
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-^4;4 

) 





^ ^o,o(i;f) , = (00000) 
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7r3/2 q 



1 1 
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1 1 - .2. 1 
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+ ^ ~7 ^6,0(1; g 



5^ 
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TT 



3/2 



+ 



^3/2 



z Ao(l;r) diag (1,1, 1,1, 1,1, 1,1,1) 

/ 14 
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-11 

70 
5v^ 
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11^ 



+ 
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which can be made partially block-diagonalized by 





( ^0,0 
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with 



'S'0,0 



'22 



'S'44 — 



V 



2^6 



1_ 




2V6 
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2V2 



]_ 

V2 





1 

V2 





1 
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1 





^\ 
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^ 
10 0/ 








VL 
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_V5_ 
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VL 
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Vl_ 
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1 
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1 

V2 
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1 

2V2 




VL 
2V2 






After this partial-diagonahzation, finite-volume function becomes F 
where 



2V6 * 

1 




2^6 




J 

(FV) 



(FV) 
0,0 

0;4 
2:2 



+ 



7r3/2 

1 
1 
1 



q 



0:2 



(00000) 



2V^ 

Zo,o(l;a diag (1,1, 1,1,1) 



- Z4,o(i;f) ^(000001000) 
-^o.ofl; "^^ 

4 ^4,o(l;g') ^ diag (-2, -2, -2, 3, 3) 



(All) 



2;4 



1 



7J-3/2 



24,o(i;r) 



20^3 
77 



;;2'\ 
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40v/3 
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^ 4;4 



1 1 



3 Ao(l;g') diag (1,1, 1,1, 1,1, 1,1,1) 



7]-3/2 q 

+ ^TTT 4 Ao(l;g') tI4 diag (2, 2, 7, 7, 7, 14, -13, -13, -13) 



7]-3/2 g5 



1001 



+ 47:7 4 ^6,o(l;g') 7-4^^ diag (-16, -16, -1,-1, -1,20, 5, 5, 5) 



7]-3/2 g7 



+ 



1 



7r3/2 g9 



^ ^8,o(i;g 



11V13 
392 



diag (7, 7, -8, -8, -8, 10, 0, 0, 0) 



(A13) 



1001V17 

It is clear from the form of the matrix F^^'' that the ordering of the F*-*) in the F4^^''-block 
is , T^, Af and T2, respectively, and the equations that dictate the energy-eigenvalues of 
each of the given in eqs. (tel), (20), (25), and (29), follow directly from these expressions. 



Appendix B: Z/ m(l;Q^) Functions 

The two-hadron Green functions in the finite lattice volume depend upon summations over 
plane-wave states subject to periodic boundary conditions and with amplitudes that depend 
upon the strength of the interactions in each of the partial-waves that generate the two- 
hadron T- matrix. The summations that define the energy-eigenvalues in the volume are [H E] 



n 



g2 



(Bl) 



a special case of the sums defined in eq. ([s]). The I = summation is special as it requires UV 
regulation in order to be defined, while sums with / > 1 are finite due to contribution from 
the solitary Yim- However, brute-force evaluation of the sums is quite inefficient and Liischer 
presented a method to evaluate the sums that exponentially accelerates their evaluation [H 
[2], making use of the Poisson resummation formula. In this appendix, we reproduce Liischer's 
results, and then present each of the Zi^m(l;g^) that contribute to the energy-eigenvalues 
considered in the body of this paper. 

Numerical evaluation of the ^o,o(l; 9^) can be evaluated by brute force through the defi- 
nition 



-= lim 

/An An-!>00 



An -I 

^ np - 



g2 



47rA, 



(B2) 
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or through the exponentially-accelerated relation [H [21 [6T] 



/\/7r ^ — ^ 



-|np 



|n|2 — g2 



E 



e * 



For / 7^ 0, the exponentially accelerated evaluation can be accomplished with 



-A 

There are exact relations that exist between the ^/ „j(l; g^) for fixed /: 



^4,±4(i; g^) 
^6,±4(l; g^) 




26,o(i;r) 



^8,±4(i;g) = — — ^8,o(i;g) , ^8,±8(i;g) = ^8,o(i;g 



33 



66 




^io,o(i;9 ) , 2:io,±8(i;9' 



;;2\ 



^i2,o(i;g 



;;2^ 



42 



-2io,±4(i; g^) = 
-2i2,±8(i;g^) = 

^i2,±i2(l;g') = 4W-^ Zi2,o(l;g') + 9J-^ Z^^Mhq 



187 

130 



^10,0 



429 
646 



4W^^i2,±4(l;g') 



91 - .2^ , n / 11 ^ n..2^ 



7429 



7429 



The Poisson resummation formula 



Y: S^y-n) = Y: e'^— - 



has been used in obtaining eq. ( B4 1 . 

We have used a relation that is similar to that used by Liischer jl], 



7r\3/2 ZTT 
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^14,±4 
-2l4,±12 

^16,±8 
^16,±12 
2^16,±16 

2^18,±8 
>2i8,±12 
^18,±16 



1;^ 
1;^ 



3 /143 



2 V 595 



214,0(1; g ) , 2i4,±8(l;g 



741 
1190 



2i4,o(i;g ) 



1 /437 



2 V 119 



2i4,o(i;g ) 



6 



442 
2185 



^i6,o(l;g ) 



31 / 13 ^ 16 / 17 ^ 



4wJ^Zi6,4(i;g') + 7 



74865 



493 
135470 



;:2\ 



l;q 



x2\ 



58 / 22 
161 



^18,4(1; g 



3 /646 



5V 23 



2i8,o(i;g 



;;2^ 



;;2\ 



l-q 



501 / 11 ^ 16 /323 ^ ,^ 



^ , 3162 ^ ,^ 



19437 
6670 



^18,0 (1;<? 



;;2\ 



(B7) 



Unlike the cases of / = 0, 4, 6, 8, 10, 14 which have only one occurrence of the irrep in 
their decomposition, / = 12, 16, 18 have two, and as such the 2^i2,±4fc {k is an integer) are 
not simply proportional to ^12,0, as demonstrated in eq. (B7), and a similar statement can 
be made about Zi6,±4fc and Zi8,±4fe. 

In an effort to better understand the origins of the structure of the functions determin- 
ing the energy-eigenvalues of each of the r*^*\ it is useful to explicitly display the functions 
Z/„i(l;g^)- The function Zo,o(l;g^) is shown in the body of this paper in fig. [ij As dis- 
cussed by Liischer, the functions 2;_m(l; g^) vanish for all odd-Z, and also vanishes for / = 2. 
The function ^4^o(l; g^) is shown in fig. ^6| and exhibits some structure that is not present 
in 2o,o(l;g^)- There are branches of 24.o(l;g^) that are non-monotonic, for instance, be- 
tween g^ = 1 and (p = 2. This behavior is found in all of the ^/.m's with I > 0. The 
functions Zs^^l^q"^) and Zioo(l;g^) are shown in fig. 

fig.g ^i4,o(i;g') in fig.|39l Zi6,o(i;g') and Zi6,4(i;g'! 

and Zi8,4(l;g^) 



in fig. 



41 



37 



■'1,'m 

2i2,o(i;g') 



and Zi2,4(l;g^) in 
in fig. 40[ and finally Zi8,o(l;g^) 



In constructing the perturbative expressions for the energy-eigenvalues in terms of the 
61, the leading contributions result from the residue of the pole of the leading function. We 
present a few of these residues of the Zi^mi^', g^) functions in table XXVIII and table XXIX 



Appendix C: Perturbative Expressions 

In many instances the energy-shifts due to the interactions are small because the phase-shift 
is small and/or the lattice volume is large. In such instances, a perturbative expression 
can be used to extract the phase-shift from an energy-eigenvalue instead of solving the full 
expression, as discussed by Liischer [HE]. The energy-eigenvalues for a given F^*) in a given 
|np-shell, and more specifically g^, can be expanded in terms of the dimensionless quantity 
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2 4 6 8 10 12 14 



FIG. 36: The function ^4^0(1; 9^) (top-panel) and ^6,o(l;^^) (bottom-panel). 

case of a single partial- wave, the general form for the energy 

^ The corresponding expansion that is appropriate for systems near unitarity, an expansion in terms of 
{LqcoiS) for s-wave interactions that is small for large scattering lengths, can be found in Ref. [56] . 
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FIG. 37: The function Z8^o(l;y ) (top-panel) and Zio,o(l;? ) (bottom-panel). 



of an irrep, r^*^ in the |n|^-shell is 



27r\ 



21 



a 



(|np,rW) 



5q 



n2 



, (CI) 
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I— I 
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FIG. 38: The function Zi2,o(l;y ) (top-panel) and Zi2,4(l;g ) (bottom-panel). 
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FIG. 39: The function Zi4,o(l; r)- 



TABLE XXVIII: Residues of the functions \/47r Zi.m{l; q^) for / < 12, i.e. the coefficient of 
where = |np + Q^. 



n ■ 



R yjA-KZ, 



--0,0 



R [^AttZqa 



R [V47r^8,oJ R [V47rZio,o 



R V47rZi2n 



R V47r^i2,4 




1 
2 
3 
4 
5 
6 

14 



1 

6 

12 

8 

6 

24 
24 
48 





21 
2 

-21 

-84 

168 

210 
-378 
-4116 





3v^ 
4 

_3!h/l3 
2 

48\/l3 
48^/l3 

-255VI3 
333\/T3 

-1806\/13 





32 
891^17 
16 

99\/l7 

792^/17 

13365v^ 



51381V17 



455301V17 





65v^ 
64 

65v^ 
32 

-520\/2l 
1040\/2T 

204425^21 
16 

408915^21 
16 

19718335^21 





3715 

256 
43885 

128 
2225 
2 

59440 

21172895 
64 

58441035 
64 

2472346835 
32 





75V1001 

512^ 

8085V1001 

256 

555V1001 
4 

600V1001 

1262535V1001 
128 

324765V1001 
128 

309589035^1001 



64 



where the solutions to eq. (CI) can be written as 

2 



° [(Lg)2'+icot(5/]' ' 



tan^5f + .. 



(C2) 
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FIG. 40: The functions Zi6,o(l;9 ) (top-panel) and Zi6,4(l;g ) (bottom-panel). 



where the phase-shift is evaluated at the unperturbed energy of the state, and the coefficients 



are 



^(|np,rW) 



^(|np,rW) 



a 



(|np,rW) 



27r2|n|2'+i 



(|n|2,rW) 



a 



(|n|2,rW) 



27r2|nP'+i 



(|n|2,r«)x2 , ^,(|n|2,r«) (|n|2,rW) 



27r2|n|2'+i 



(C3) 
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FIG. 41: The functions Zisfl{l;q ) (top-panel) and ^18^4(1; g ) (bottom-panel) 



For the terms in eq. (C2), the contributions scale as 



1^61+3 



O 



(C4) 



from which it can be determined when the contributions from higher partial-waves become 
important. For instance, the energy-shifts in the T{ irrep are dominated by 5i, and the 
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TABLE XXIX: Residues of the functions V47r Zi^^il-^q^) for 14 < / < 18, i.e. the coefficient of 
where q'^ = |np + Q^. 



|n|2 


R \/^Zi4fl 


R [VA^Ziefi 


R [V4^Zi6,4] 


R [\/47rZi8,o] 


R -v/47r2^i8,4 


1 

2 
3 
4 
5 
6 

14 


595^ 


22819^33 


33^323323 


20613v^ 


39^920227 


512 
52955v^ 
256 

3570\/29 
19040\/29 

9356375^29 
128 

144911655v^ 


8192 
2627491733 
4096 
135261 


4096 


16384 
234021v^ 
8192 

-12453^/37 
329808\/37 

2195267925^37 

4096 
40787860977 


8192 


5937^323323 

2048 
873V323323 


9063^920227 
4096 

18^920227 


16 

182552\/33 

3105175795^33 

2048 
1849263939^33 


8 


528\/323323 


-1248\/920227 


1267665^323323 
1024 


12933225^920227 
2048 


15619473^323323 


26504469V920227 


128 

6869655205^29 


2048 

12080581399901^33 


1024 


4096 

44894607368667^37 


2048 


4554268593V323323 


51258039399V920227 


64 


1024 


512 


2048 


1024 



(InP A'^) 

TABLE XXX: The coefficients, ' ^ that contribute to the perturbative expansion of the 
energy-eigenvalues of states in the Af irrep of the cubic group, as given in eq. (C2), in terms of 
s-wave phase-shift Sq. 



w 

1 

2 
3 
4 
5 
6 



% 

3 

n'^ 
4 

V37r2 

3_ 

27r2 
12 

_2x/6 

^2 



fh 

-0.06137 

-0.1826 

-0.1981 

0.2415 

0.1590 

-0.4798 



^2 

-0.3542 
-0.3618 
-0.1996 
-0.1328 
-0.5155 
-0.2025 



expansion is of the form 

~ Oil) ^o{L)^ + O(^) + O(^) + ... . (C5) 

respectively, and the I = 3 partial- wave first contributes at O (jy)- In the case of the 
irrep, which is dominated by ^4, the expansion is of the form 

^„P,T. ~ 0(1) + '^{t') ^ "^(is) ^""{t^) ^ ""{t^) + - • (°«' 

and ^6 contributions are of the form 0{L~^^). Therefore, the order at which the higher 
partial- waves contribute in the large- volume limit depends upon the T^^\ 

The perturbative expansions of the lowest few energy-eigenvalues in terms of the / = 
phase-shift 60 were given by Liischer pQ, and here we simply extend those results to levels 
with |np < 6 with the coefficients given in table 

|np = level can be expressed in terms of the s-wave scattering parameters defining the 



XXX The energy of the Af state in the 
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(Inp T"'"""') 

TABLE XXXI: The coefficients, that contribute to the perturbative expansion of the 

energy-eigenvalues of states in the (first occurrence of the) Tj~ irrep of the cubic group, as given in 
eq. (C2), in terms of / = 1 phase-shift 6i. 



\n\' 
1 
2 
3 
4 
5 
6 



% 

__3_ 
_3V2 



3 
2^ 
12 
V57r2 
205 



!h 

-0.3653 

-0.3975 

-0.2761 

0.2035 

0.05024 

-0.5625 



92 



-0.2058 
-0.1979 
-0.1471 
-0.1589 
-0.5555 
-0.07659 



TABLE XXXIL The coefficients, g}' ' ' ' that contribute to the perturbative expansion of the 
energy-eigenvalues of states in the (first occurrence of the) irrep of the cubic group, as given in 
eq. ( C2 ) , in terms of / = 2 phase-shift 82 ■ 



|n|2 




(|n|2,£;+(i)) 

9i 


(|n|2,i?+(i)) 
92 


1 


15 
2^ 


-1.5672 


2.5842 


2 


15 
4v^7r2 


-0.8065 


0.54 


4 


15 
47r2 


0.3272 


-0.421 


5 


78 
5v^7r2 


-0.447 


-0.4331 


6 


5 /s 
2^\J 2 


-0.7884 


0.3746 



low-energy behavior of the phase-shift, and it is well-known that 



^2 



+ ... 



(C7) 



where the particle-physics convention for defining the scattering length has been used, and 
the coefficients are Ci = —2.8373 and C2 = 6.3752. 

An impo rtant point to note is that the perturbative energy-shifts that are presented in 
table XXX-table XXXIII are for one of the occurrences of the F*-*) that form a given |np- 
shell. Other occurrences are unperturbed at leading order. When multiple occurrences of 
a given irrep appear in a given |np-shell, the leading interactions will perturb the energy 
of one combination, while leaving the other states unperturbed, but the interactions in 
higher partial-waves will perturb these remaining states. The expansion coefficients for the 
lowest-lying T{ (dominated by the and T2 (both dominated by 5|) are shown in 



table XXXI table |XXXTT| and table XXXIII respectively. We note that the coeffici ents in 
the perturbative expansion of the energy-eigenstates in the irrep given in table XXXI 
differ from those given by Liischer [2]. This can be attributed to the fact that g^^^+^-'cot^; 
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TABLE XXXIII: The coefficients, 5-'"'''^"^ 



(1) 



that contribute to the perturbative expansion of the 
energy-eigenvalues of states in the (first occurrence of the) irrep of the cubic group, as given in 
eq. (C2), in terms of I = 2 phase-shift 62- 



\n\' 
2 
3 
5 



9o 



(|np,T+('') 



15 



2v^7r2 

20 
3v^7r2 

48 

°V 2 









92 


-0.4830 


-0.1828 


-0.6737 


0.2128 


0.2004 


-0.4515 


-0.5497 


-0.0902 



(|n|2 r) 

TABLE XXXIV: The coefficients, Qn that contribute to the perturbative expansion of the 



energy-eigenvalues of states in the A2 , T2 , , E~ and as given in eq. (C2), in terms of 
dominant phase-shifts iJs, (Js, ^4, iJs, and Je, respectively. Also given are the coefficients in the 
perturbative expansion of the second occurrence of E 
phase-shifts 5^, ^4, 5^, and ^5, respectively. 



rp+ rp- 



and T2 , in terms of dominant 



|n|2 


9o 


9o 


(|n|^T'+) 

yo 


9o 


9o 


(|n|2,B+(2)) 

yo 


(|n|2,T+'^)) 

yo 


(|n|2,T-(2') 

yo 


(|n|2,T-(2') 

yo 


2 
3 
5 

6 




140 


105 





2268 






385 





162162 





23814 






245v^ 





252 





74844 


8v^7r2 



84 





5V57r2 
35 


125V57r2 
35v^ 


3125V57r2 




1625V57r2 




25V57r2 
175v^ 


3125V57r2 




3tt2 


4v^7r2 


4v^7r2 


12v^7r2 


36v^7r2 


36v^7r2 



can be expanded in a power-series in energy about threshold as performed in this work, 
while qcotSi does not have such an expansion for / > 0. For the remaining irreps, the 
T2 , A2 , T^, A2 and E~ , the expansion converges rapidly with just one non-trivial term, 
O (tan . The leading coefficients for these expansion of the energy-eigenvalues for each 
of these irreps are given in table XXXIV, along with the coefficients in the expansions for 

Tj_and Tg". The perturbative expansion of the 



the second occurrences of the i?"*", T/, 



lowest-lying state is given in table XXVII 
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